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Abstract. IRAS 08339+6517 is a luminous infrared and Lya-emitting starburst galaxy that possesses a dwarf 
companion object at a projected distance of 56 kpc. An HI tidal tail has recently been detected between both 
galaxies, suggesting that about 70% of the neutral gas has been ejected from them. We present deep broad-band 
optical images, together with narrow band Ha CCD images, and optical intermediate-resolution spectroscopy of 
both galaxies. The images reveal interaction features between both systems and strong Ha emission in the inner 
part of IRAS 08339+6517. The chemical composition of the ionized gas of the galaxies is rather similar. The 
analysis of their kinematics also indicates interaction features and reveals an object that could be a candidate 
tidal dwarf galaxy or a remnant of an earlier merger. Our data suggest that the Hi tail has been mainly formed 
from material stripped from the main galaxy. We find weak spectral features that could be attributed to the 
presence of Wolf-Rayet stars in this starburst galaxy and estimate an age of the most recent burst of around 
4-6 Myr. A more evolved underlying stellar population, with a minimal age between 100 - 200 Myr, is also 
detected and fits an exponential intensity profile. A model which combines 85% young and 15% old populations 
can explain both the spectral energy distribution and the Hi Balmer and Hel absorption lines presented in our 
spectrum. The star formation rate of the galaxy is consistently derived using several calibrations, giving a value 
of ~9.5 Mg yr -1 . IRAS 08339+6517 does satisfy the criteria of a luminous compact blue galaxy, rare objects 
in the local universe but common at high redshifts, being a very interesting target for detailed studies of galaxy 
evolution and formation. 


1. Introduction 

Determining the star formation history in galaxies is fun¬ 
damental for the understanding of their evolution. In 
nearby galaxies, such as those belonging to the Local 
Group, it is performed studying the stellar content via 
color-magnitude diagr ams (CMD) and stars as age trac¬ 
ers [see ICrebell (1999) and references there in]. However, 
in distant galaxies it is only interpreted by means of spec¬ 
tral synthesis techniques. The analysis of the interstellar 
medium (ISM) in and around galaxies complements the 
study of their star formation histories. Star formation ac¬ 
tivity is stronger in starburst galaxies, where the domi¬ 
nant and very young stellar population could even hide 
the more evolved stellar population presented in it. The 
more extreme cases are the blue compact dwarf (BCDs) 
galaxies, that have very low metal abundances and ex¬ 
hibit a global starburst activity. Although the majority 
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* Based on observations made with NOT (Nordic Optical 
Telescope) operated on the island of La Palma by Instituto de 
Astroffsica de Canarias in the Spanish Observatories of Roque 
de Los Muchachos of the Instituto de Astroffsica de Canarias. 


of these systems possesses an o ld u nder lyi ng pop ulation 
with ages of several Gyrs llLoose fc Thuanifl986 L a few 
of them do not p resent it, suggesting t hat they are really 
young galaxies lllzotov fc Thuaul Il999i l. The best exam¬ 
ple of this kind of BCDs is I Zw 18, a bona fide young 
galaxy in the local universe lllzotov fc Thuaul 120 041. At 
intermediate and high redshifts a heterogeneous class of 
vigorous starburst systems with lu minosities arou nd L* 
( L* = 1.0 x 10 1o Lq) are observed ( Guzman et al ] 2001. 
They are designated as luminous compact blue galaxies 
(LCBGs) and their evolution and nature are still open 
questions, being fundamental to get a sample of them in 
the local uni verse to investigate the origin of their activity 

( Werk et IDEnoi. 


An important subset of starbursts are the so-called 
Wolf-Rayet (WR) galaxies, whose integrated spectra have 
a broad emission feature around 4650 A that has been at¬ 
tributed to WR stars. This feature consists of a blend 
of N ill A4 640, C III/ C iv A4650 and He n A4686 emis¬ 
sion lines (Conti: 19911. the last one being the most 
prominent line. The WR feature indicates the pres¬ 
ence of a substantial population of these kinds of mas¬ 
sive stars whose ages are less than 6 Myr and of- 
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fers the opportunity to study very young starbursts 
dSchaerer fc Vaccal 1 1 998ti . Furthermore, they constitute 
the best direct measure of the upper end of the initial 
mass function (IMF), a fundamental ingredient fo r studw 
ing unresolved stel lar populations llSchaerer et al.l l2000: 
IPindao et al.l l2002tl. Stu dying a sample of WR galax- 
ies. iMendez fc Estebanl 1200(11 suggested that interac¬ 
tions with or between dwarf objects could be the main 
star formation triggering mechanism in dwarf galaxies 
and noted that the interacting and/or merging nature of 
WR galaxies can be detected only when deep and high- 
resolution images and spectra are a vailable. The compila- 
tion of WR galaxies performed bv ISchaerer et al.l ( 1999 1 


creased (PoDescu & HoduI 200(i 

Berevall & Ostlinl 

2002: 

Contini et al. 20od Pindao et alJ 

2002 

Lfpari et al. 

2003: 

Tran et al. 2003; Fernandes et al. 2004 

Izotov et al. 2004: 

Pustilnik et al. 20041 Kniazev et al. 20041). these galaxies 


12004+ In this paper, we add a new member to the list of 
probable WR galaxies: the luminous infrared galaxy IRAS 
08339+6517, which, furthermore, could be also classified 
as a luminous compact blue galaxy (LCBG). 

IRAS 08339+6517 is at 80 Mpc (at that distance, 
1" = 388 pc, assuming Ho=75 km s -1 Mpc -1 ) and was 
first ly reporte d in the IRAS Point Source Catalog (1986). 
iMargon et ali lj 19881 1 performed a multi-wavelength study 
and described the object as an exceptionally bright 
and compact starburst nucleus. iGonzalez-Delgado et alJ 
ljl998tl presented a detailed UV analysis, suggesting that 
O stars must be present in the ionizing cluster(s) of the 
galaxy in order to explain the blue wing of the Civ 
A1550 and SilV A1400 stellar absorption lines, although 
these stars cannot be more massive than about 40 M 0 . 
iKunth et al.l ( 1998 ! noted a P Cygni profile in the Lya 
emission, which also showed two components indicating 
the c haotic structure of the i nterstellar medium in this ob¬ 
ject. iMas-Hesse et, al.l (l2003ll included IRAS 08339+6517 
in their HST UV and optical stud y of Lya s t arbur sts. 

In a very recent paper, ICannon et al. ll2004ll pre¬ 
sented VLA Hi imaging of this starburst galaxy and 
found an extended tidal structure in neutral hydrogen 
indicating that it is interacting with a nearby compan¬ 
ion, 2MASX J08380769+6508579 (which we designate as 
the companion galaxy throughout this paper). This fea¬ 
ture provi des an ev idence for tidally induced starburst 
episodes. ICannon et al. 1 12 0041 ) estimated neutral hydro¬ 
gen masses of (1.1 ± 0.2) x 10 9 Mq and (7.0 ± 0.9) x 10 s 
Mq for IRAS 08339+6517 and its companion galaxy, re¬ 
spectively, and a mass of (3.8 ± 0.5) x 10° Mq in the tidal 
material between them. Consequently, ~ 70% of the neu¬ 
tral gas has been removed from one or both galaxies. 

We are developing a detailed analysis of Wolf-Rayet 
galaxies using both deep images and spectra in optical 
wavelenghts. Our data show increasing evidences of the 
connection between the detection of WR features and the 
findi ng of interaction s ignatu re s in this kind of starbursts 
l|Mendez fc Esteban! 2000 : iLonez-Sanchez & Estebanl 


Table 1 . Summary of observations performed using 
ALFOSC in 2.56m NOT, with a spatial resolution of 0.19" 
pix -1 . The position angle, spectral resolution and wave¬ 
length range of the intermediate resolution spectroscopy 
is shown in the last line. 


Observations 

Date 

Exp. Time (s) 

Filter 

Broad-band 

05/04/03 

3 x 300 

U 

imaging 

05/04/03 

3 x 300 

B 


04/03/20 

2 x 300 

V 


04/03/20 

3 x 300 

R 

Narrow-band 

04/03/20 

3 x 300 

Ha (F.19) 

imaging 

04/03/20 

2 x 300 

Ha cont. (F.20) 

Intermediate 

04/03/20 

3 x 900 

grism #7 

resolution 

P.A. 

R 

AA 

spectroscopy 

138° 

1.5 A pix -1 

3600-6800 A 


2003; iLonez Sanchez et al l2004alht Lopez-Sanchez 
2006+ The characteristics previously observed in IRAS 
08339+6517 (strong Ha emission, starburst nature and 
the detection of an extended H 1 tidal tail in the direction 
of its companion dwarf galaxy) make it an ideal target to 
look for the existence of WR stars. In this way, we have 
performed deep broad-band optical and narrow-band 
Ha CCD images together optical intermediate-resolution 
spectroscopy of IRAS 08339+6517 and its dwarf compan¬ 
ion galaxy. Our main aim is to study their morphology 
and stellar populations, as well as the distribution, 
chemical composition and kinematics of the ionized gas, 
and to check if WR stars are presented in their youngest 
bursts. 


2. Observations and data reduction 

The main data were obtained on 2004 March 20 at the 
2.56 m Nordic Optical Telescope (NOT) at Roque de 
los Muchachos Observatory (La Palma, Canary Islands, 
Spain). We used the ALFOSC (Andalucfa Faint Object 
Spectrograph and Camera) instrument in both image and 
spectrograph mode, with a Loral/Lesser CCD detector 
(2048 x 2048 pixels) with a pixel size of 13.5 /im and 
spatial resolution of 0.19" pixel -1 . Additional U and B 
images were obtained on 2005 April 03 using the same 
telescope and instrument. The journal of all the imaging 
and spectroscopical observations can be found in Tabled 


2.1. Optical imaging 

We obtained three 300 second exposures for U, B and 
R standard Johnson filters and two 300 second exposures 
for the V filter that were combined to obtain the final 
images with a P SF FWH M of 0.8". The standard starfield 
SA107-602 iLandohlll 99l was used to flux calibrate the 
final images. It was observed at different airmasses during 
the night. Twilight images of different zones of the sky 
(blank fields) were taken for each filter in order to perform 
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Fig. 1 . (a) Deep optical image of IRAS 08339+6517 and its companion galaxy, 2MASX J08380769+6508579, in R. It has been 
saturated to reveal the weakest features, (b) The companion galaxy and (c) IRAS 08339+6517 in non-satured R images. The 
grayscale is in logarithmic scale in all three cases. 



the flat-field correction. All the reduction process (bias 
subtraction, flat-fielding and flux calibration) was done 
with IRAF 1 package following standard procedures. IRAF 
software was also used to determine the 3a contours over 
the sky background level and to obtain the photometric 
values of the galaxies. 

2.2. Ha imaging 

The narrow-band filters for Ha and red-continuum ima¬ 
ges (centered at 6687 and 6571 A and with a FWHM of 
50 and 47 A, respectively) were selected taking into ac¬ 
count the recession velocity of the object given by the 
NASA/IPAC Extragalactic Database (NED). Three (two 
for the continuum) exposures of 300 seconds were added to 
obtain the final image in each filter. The reduction process 
was similar to the broad-band images. The absolute flux 
calibration was achieved by taking short ex posures of the 
spect rophotometric standard star Feige 56 dl lamuv et al.l 
11992) at different airmasses. 

To obtain the continuum-subtracted Ha image, we fol¬ 
lowed the standard procedure given bv [Cedresl ( 200 3). A 
scaling factor between the Ha and the continuum frames 
was determined using non-saturated field stars (in our 
case, this scaling factor was 0.9). The FWHM of the 

1 IRAF is distributed by NO AO which is operated by AURA 
Inc., under cooperative agreement with NSF 


PSF for both images was similar (around 0.6"), so no 
Gaussian filter was applied. The finally flux-calibrated 
and continuum-subtracted Ha image was then produced 
by subtracting the scaled continuum frame from the Ha 
frame. The W(Ha) image was also obtained following the 
standard procedure, correcting the W(Ha) for the under¬ 
lying stellar population of the disk. 


2.3. Intermediate-resolution spectroscopy 

We used ALFOSC in spectrographic mode with a slit of 
6.4' (the field of view of the camera) long and 1" wide to 
obtain intermediate-resolution spectroscopy of the galax¬ 
ies. Grism #7 with 600 rules mm -1 and resolution of 1300 
was used, with a dispersion of 111 A mm -1 and spectral 
resolution of 1.5 A pix -1 . Each spectrum spans the wave¬ 
length range 3200 to 6800 A (redsliift corrected). The slit 
position was located at P.A. = 138°, which was chosen in 
order to cross the center of IRAS 08339+6517 and its com¬ 
panion galaxy. Three 15 minute exposures were taken and 
combined to obtain good signal-to-noise and an appropri¬ 
ate removal of cosmic rays in the final spectra. We used 
comparison He and Ne lamp exposures for the wavelength 
calibration of the spectra. The absolute flux calibration 
was ac h ieved by observ ations of the standard star Feige 
56 dMassev et alill988 ). The observations were made at 
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Table 2. General properties of the galaxies 


Property 

IRAS 08339+6517 

companion 

M b ° 

— 21.57 ± 0.04 

-18.21 ± 0.04 

ms° 

12.94 ± 0.04 

16.31 ± 0.04 

L b/Lq 

(6.69 ± 0.24) x 10 10 

(2.99 ± 0.11) x 10 9 

E (B - V+ 

0.16 ± 0.02 

0.13 ± 0.02 

(U - B)° 

-0.51 ± 0.08 

-0.16 ± 0.10 

(B - V)° 

0.02 ± 0.08 

0.20 ± 0.08 

(V - R)° 

0.24 ± 0.08 

0.26 ± 0.08 

(V - J)°’ h 

1.39 ± 0.06 

1.56 ± 0.12 

(J — H) 0,b 

0.64 ± 0.05 

0.21 ± 0.25 

(.H - K s ) 0,h 

0.23 ± 0.06 

0.68 ± 0.28 

Mhii (Mq) 

(1.45 ± 0.07) Xl0 7 

(2.6 ± 0.3) x 10 5 

M H i (Mq) C 

(1.1 ± 0.2)xl0 9 

(7.0 ± 0.9) X10 8 

Mdust (Mg) d 

4.5 Xl0 6 


M Kep (Mq) 

(10 ± 3) X10 9 

(8 ± 2) x 10 9 

MfCep/Mni 

9.1 

1.1 

Mrii/Mhi 

0.013 

0.0004 

Lb/M Kep 

6.7 

0.37 

Lb/Mr i 

60.8 

4.2 

L B /M dust 

1.5 xlO 4 


SFR (Mq yr" 1 ) 

9.5 

0.17 

Tgas (Gyr) 

0.15 

5.4 

Age e (Myr) 

4-5 

5-6 

Age f (Myr) 

30 - 50 

> 250 

Av r g (km s -1 ) 

0 

20 ± 10 

12+log O/H* 1 

8.45 ± 0.10 

8.38 ± 0.10 

N(WR)‘ 

~ 310 


WR/(WR+0)‘ 

0.03 



a Derived from C(H (3). See §3.1 for more details. 
b Derived from the values of 2MASS shown by the NED. 
c From lCannon et al. I I 2004)_ 

d Derived using the IHuchtmeier et al.l J1995I) relation between the 60 
and 100 /im FIR fluxes and the dust mass (last equation in p. 678). 
e Estimated age of the most recent bursts using our Ha and 
spectroscopical data. See §4.1 for more details. 

f E stimated age compari ng our optical col ors with the STARBURSTS 
99_ JLeitherer *eT HJEEaH) and PEGASE.2 {F ioc Sz Rocca-Volmerangd 
I1997T models. See §4.1 for more details. 

s Radial velocity wi th respec t to the center of IRAS 08339+6517. 
h Derived using the lPilvuginl J200U) empirical calibration (see §3.4.2). 

1 Approximate number of WR stars and WR/(WR+0) ratio (see §4.3). 


air masses between 1.1 and 1.3. Consequently, no correc¬ 
tion was made for atmospheric differential refraction. 

All the CCD frames for spectroscopy were reduced us¬ 
ing standard IRAF procedures to perform bias correction, 
flat-fielding, cosmic-ray rejection, wavelength and flux cal¬ 
ibration, and sky subtraction. The correction for atmo¬ 
spheric extinction was performed using an average curve 
for the continuous atmospheric extinction at Roque de los 
Muchachos Observatory. For the two-dimensional spectra 
two different apertures were defined along the spatial di¬ 
rection to extract the final one-dinrensional spectra of the 
galaxies, centering each aperture on its brightest zone. 
Small two-dimensional distortions were corrected by fit¬ 
ting the maxima of [Oil] AA3726,3729 doublet and Ha 
emission lines (between 2 and 3 pixels, 0.4 - 0.6"). 


3. Results 

3.1. The extinction toward IRAS 08339+6517 

We have carefully studied the reddening coefficient, 
C(H/3), derived from our spectra using the Balmer decre¬ 
ment (see §3.4) in order to achieve accurate photomet¬ 
ric and spectroscopic data of IRAS 08339+6517 and its 
companion galaxy. However, the fluxes of nebular Balmer 
lines in the optical spectra must be also corrected for un¬ 
derlying stellar absorption. In this way, we performed an 
iterative procedure to derive the reddening coefficient and 
the equivalent widths of the absorption in the hydrogen 
lines, W a 6 s , to correct the observed line intensities for both 
effects. We assumed that the equivalent width of the ab¬ 
sorption components is the same f or all the Balmer lines 
and used the relation given by iMazzarrella fc Boronsonl 
(Il993i) to the absorption correction, 


c{m = W )'° s 


RA) 

nm 


X 1 + 


W a , 


Wh/3 J 


F+) 

F{Hp) 


( 1 + 


W„, hs 

w x 


( 1 ) 


for each detected hydrogen Balmer line. In this equation, 
F(A) and /(A) are the observed and the theoretical fluxes 
(unaffected by reddening or absorption); W a bs , Wa and 
Wh/3 are the equivalent widths of the underlying stel¬ 
lar absorption, the considered Balmer line and H j3, re¬ 
spectively, and /(A) is the reddening curve normalized to 
H 8 using the ISeaton Il97f+ law. We have used the pairs 
Ha/H/3, Hy/H/J and Hd/H/3 for the main galaxy but only 
the two first line ratios for the companion object. We have 
considered the theoretical ratios between these pairs of H 1 
Bal mer lines expec t ed for case B recombination given by 
IStorev fc Hummeii l|l995l) and appropriate electron tem¬ 
peratures around 10 4 K and electron densities around 100 
cm" 3 (Ha/H/3 = 2.86, H+/H/3 = 0.469 and HS/U0 = 
0.259). The C(H/3) and W a b s that provide the best match 
between the corrected and the theoretical line ratios are 
shown in Table fl] 

Once we have obtained the reddening coefficient, the 
photometric data can also be corrected for ext inctio n by 
interstellar dust. We used the standard ISeatonl ( 1979 ) law 
and the standard ratio of Ry = Ay/E(B — V) = 3.1 to 
obtain the extinction in V, Ay ~ 2.11 x C(H/3) and the 
color excess, E(B — V) = 0.68 x C(H/3). The data in U , 
B and R were corrected using Ay = 1 531 x Ay , A r = 
1.324 x Ay and Ar = 0.748 x Ay Jliieke fc Lebofskvl 
Eiil, whereas the data in Ha filter were corrected using 
Ah a. = 1.61 x C(H/3) llSeatonlll 979li . 

As we explain in §3.4, we have extracted eight different 
apertures to study the ionized gas in IRAS 08339+6517. 
The C(H/3) derived varies between 0.18 and 0.30, imply¬ 
ing a color excess between 0.12 and 0.20. We have adopted 
the average C(H/3) (= 0.24) of these eight zones as rep¬ 
resentative for the galaxy. The derived color excess used 
to correct the photometric data is E(B — V) = 0.16 (see 
TablcQ). 
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Fig. 2. Continuum-subtracted Ha image of IRAS 
08339+6517. The contours represent 3(7, 5a, 7<r, 15cr, 
30cr, 70 ct, 140a, 230(7, 430(7 and 700cr above the sky level. We 
also show the position of the slit used for spectroscopy, with 
P.A. 138°, and label the different individual regions that we 
have analyzed. 



Fig. 3. Relative Ha flux along the slit (see Figure^. The dif¬ 
ferent zones for which we have extracted spectra are indicated 
and labeled. Knots #1 (the center of the galaxy), #2, A and 
B are 1" in size. We also show four additional apertures ex¬ 
tracted to study spatial variations of the ionized gas inside the 
galaxy. They have 2.5 /x , 5", 9” and 12" in size. The aperture 
of 12" represents the integrated spectrum of the galaxy, that 
we show in Figure Pf top). 


ISchlegel et al.l ( 1998 ) presented an all-sky map of in¬ 
frared dust emission, determining a color excess E(B — V) 
= 0.092 (assuming Ry = 3.1) at the Galactic longitude 
and latitude of IRAS 08339+6517. This value is lower than 
the E{B — V") we obtain for the galaxies (0.16 and 0.13), 
indicating that a considerable fraction of the dust in the 
line of sight is external to the Milky Way, especially in the 
main object. 

The importance of t he extinction co rr ectio n was pre¬ 
viously remarked on bv IGonzalez-Delgado et all ( 1998 1. 
who studied the variation of the color excess derived 
from various methods: the Balmer decrement (that gave 
an E(B — V) = 0.52, a simil ar value that the one de¬ 
rived by iMargon et al.l Il98ij) using the Ha/H/3 ratio, 
E(B — V) = 0.55, although they do not correct for the 
underlying stellar absorption) and the change of the UV 
continuum slope [which indicates lower values for the ex¬ 
tinction, E(B — V) = 0.1 — 0.2]. They explain this dif¬ 
ference in terms of a nonhomogeneous distribution of the 
dust, although it could also be due to the variation of 
the evolutionary stage of the stellar population with the 
aperture size. The data derived from our optical spec¬ 
tra sug gest t h at E( B — V) is n ot so high a s those de¬ 
rived bv lGonzalez-Delgado et al ] ill 998ll and lMargon et ahl 
(Il988h using the same method but is similar to the value 
derived fro m the UV data. This result is supported by 
IStevens fc Strickland! (Il998tl . who studied ROSAT PS PC 
seredipitious observations of IRAS 08339+6517 and con¬ 
cluded that the spectral fitting yields an X-ray tempera¬ 
ture of kT = 0.58 keV, with Z = 0.02 Zq and a column of 
1.40 x 10 21 cm -2 , implying substantial local a b sorpti on. 
Using the dust-to-gas ratio given bv lKent et al.1 (. 1991 ). a 
color excess of 0.24 is derived, a bit higher than the one 


we have assumed from tire Bal mer decrement, but not as 
high as that derived bv IGonzalez-Delgado et all MUi). 

However, we have noted significant differences in 
C(H/3) derived for different zones and apertures inside 
IRAS 08339+6517. We briefly discuss this in §4.4. 

3.2. Optical imaging 

In Figure we show our final deep R image, which has 
been saturated to show the faintest structures. This deep 
image reveals a disturbed morphology in the outer regions 
of IR AS 08339+6517. The most relevant are a bright ray at 
the south of the galaxy and a long arc of material connect¬ 
ing the north of the galaxy with the southern bright ray. 
This arc seems to be broken into two parts to the SE of the 
galaxy. Another weaker arc goes from the west towards the 
end of the bright ray. Furthermore, a very diffuse plume 
is found in the NW of the galaxy, precisely in the direc¬ 
tion towards the companion galaxy. This plume coincides 
wi th the beginning of the tidal HI material discovered by 
ICannon et_al,_l (2J2Q4) between the two galaxies. 

The non-saturated images of the galaxies are shown in 
Figures^) and c. IRAS 08339+6517 shows a circular mor¬ 
phology. Assuming that it has a disk geometry, this fact 
would indicate that the galaxy is nearly face on. Several 
structures are found in its central areas following an ap¬ 
parently symetrical pattern from the nucleus. Although 
they are not clearly resolved, they show lower U — B val¬ 
ues and possess important Ha emission features (see next 
section), indicating that the inner area of the galaxy has 
relatively high star-formation activity. These facts could 
suggest the presence of a disk-like structure inside IRAS 
08339+6517. 
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Table 3. Ha properties of the galaxies and studied knots inside IRAS 08339+6517 


Knot 

Flux (10 -13 erg cnW 2 s -1 ) 

Luminosity (10 41 erg s 1 ) 

SFR (Mq yr" 1 ) 

Men (10 e M o ) 

— W (Ha) (A) 

Age (Myr) 

IRAS 08339+6517 

16.0 ± 0.7 

12.0 ± 0.6 

9.5 ± 0.5 

14.5 ± 0.7 

110 ± 10 

4.4 ± 0.2 

# 1 

4.3 ± 0.3 

3.3 ± 0.2 

2.62 ± 0.18 

4.1 ± 0.3 

150 ± 15 

3.3 ± 0.2 

# 2 

2.3 ± 0.2 

1.74 ± 0.14 

1.38 ± 0.11 

2.09 ± 0.17 

103 ± 9 

4.4 ± 0.2 

# 3 

0.47 ± 0.05 

0.36 ± 0.04 

0.28 ± 0.03 

0.43 ± 0.04 

129 ± 13 

3.9 ± 0.2 

# 4 

0.76 ± 0.08 

0.57 ± 0.06 

0.45 ± 0.05 

0.69 ± 0.07 

101 ± 10 

4.4 ± 0.2 

# 5 

0.24 ± 0.03 

0.18 ± 0.02 

0.14 ± 0.02 

0.22 ± 0.03 

95 ± 11 

4.5 ± 0.3 

# 6 

0.70 ± 0.07 

0.53 ± 0.05 

0.42 ± 0.04 

0.63 ± 0.06 

102 ± 13 

4.4 ± 0.3 

# 7 (NW arm) 

0.096+ 0.015 

0.073+ 0.011 

0.06 ± 0.01 

0.087 ± 0.013 

40 + 6 

5.4 ± 0.3 

Companion 

0.28 ± 0.03 

0.21 ± 0.02 

0.17 ± 0.02 

0.26 ± 0.03 

31 + 5 

5.5 ± 0.3 


The companion galaxy has a not very excentric ellipti¬ 
cal morphology, with no relevant feature inside it. It shows 
bluer colors in its external areas than in its center. 

We have performed integrated aperture photometry 
for the galaxies in the optical filters. We have determined 
the area for which we have integrated the flux using the 3er 
level isophote over the average sky level in the B image. 
We have used the same area in U , V and R, correcting 
all values for reddening as indicated in the previous sec¬ 
tion. We have determined the errors for the photometry 
by considering the FWHM of the PSF, sky level and flux 
calibration for each frame. 

The final integrated photometric values derived for 
each galaxy are listed in Table |3 We also show the NIR 
colors derived from the values given by 2MASS, that are 
also corrected by reddening. IRAS 08339+6517 has bluer 
colors than its companion object. A detailed analysis of 
the surface brightness reveals the presence of two differ¬ 
ent stellar populations inside IRAS 08339+6517, as we will 
see in §4.1. The B magnitudes derived implying an abso¬ 
lute blue magnitude of Mb ~ —21.6 for the main galaxy 
and Mb ~ —18.3 for its companion. 


3.3. Ha imaging 

The Ha map of IRAS 08339+6517 is shown in Figure [5] 
The slit position used for spectroscopy, as well as different 
regions inside the galaxy that we have also extracted and 
analyzed, are indicated and labeled. The galaxy has an al¬ 
most circular morphology in Ha emission, except for the 
arm located to its northwest (region #7 in Figure [5|. This 
feature is also in the direction of the companion galaxy. 
Several star-formation knots are found inside the disk, the 
central (region #1) being the brightest one. This indicates 
that the entire galaxy is undergoing a considerable star- 
formation activity. Regions $=4 (which seems to be con¬ 
nected to #1 and is broken into three small knots) and 
#6 (connected to #2) are elliptic and extend from the 
center to the outer regions, suggesting diffuse spiral arms. 
The NW arm also shows three knots of weak Ha emission. 
Only the brightest is found in the optical images, showing 

U — B -0.4, B-V ~ 0.1 and V — R ~ 0.2 colors. The 

companion galaxy is weaker in Ha emission but it also 
shows star-formation activity. Its morphology is irregular 


and split into several weak knots located not in its center 
but in its outer areas. 


We obtained the Ha fluxes for both galaxies and some 
regions inside IRAS 08339+6517 using this flux-calibrated 
continuum-subtracted Ha image. The boundaries of the 
zones were determined selecting each distinct emission 
peak as a separate region by eye. We have corrected the 
Ha fluxes for both reddening (see §3.1) and [N ii] contam¬ 
ination. The [N n] contamination was corrected using the 
[Nil] A6583/Ha and [Nil] A6548/Ha ratios derived from 
our spectra and the transmittance value of the Ha fil¬ 
ter for these wavelengths. In the zones for which we have 
no spectra we used the average [N n]/Ha ratios obtained. 
The errors for the Ha flux was determined considering the 
FWHM of the PSF, sky level, the flux calibration, and the 
contributions of the reddening and [N n] contamination, 
which vary between 4% (for the main galaxy) and 16% 
(for the NW arm). The final fluxes are shown in Table 0 
The Ha flux derived for IRAS 084 49+6517 is in excellen t 
agreement with that obtained by iMargon et alJ (1988), 


FH a = 1.3 x 10 12 erg s " cm 

The luminosity of the galaxies and regions was calcu¬ 
lated assuming a distance of 80 Mpc for the system. For 
comparison, the giant extragalactic H II region 30 Doradus 
in the LMC has an Ha luminosity of 1.5 x 10 40 erg s _1 
(|Kennicuttl Il984t) . similar to region #5. Assuming that 
the ionization is due to hot stars and that a single 07V 
star con tributes with a luminos ity of Ln a = 1-36 x 10 36 
erg s _1 I Schaerer fc Vaccall 1998ll to the total flux, we esti¬ 
mated a population of around ~88000 07V stars for IRAS 
08339+6517. 


,-i 


-2 


We have derived the star formation rate (SFR) and 
the mass of the ionized gas, Mhii , using the standard 
iKennicuttl lll998l) calibration. For IRAS 08339+6517 we 
find a SFRn a = 9.5 M 0 yr _1 , which is in excellent agree¬ 
ment with previous studies. We will compare the SFR 
derived from different wavelengths in §4.2. For the mass 
of the ionized gas, we have derived a value of Mhii ~ 
1.5 x 10 7 Mq. 

We have obtained the equivalent width of the Ha emis¬ 
sion, W(Ha), of the galaxies and the knots inside IRAS 
08339+6517 from the Ha image. The results are compiled 
in Tabled All the objects present very similar W(Ha), 
except the center of the galaxy (knot #1, which has the 
highest value) and the companion galaxy and the NW 
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Fig. 4. Redshift-corrected spectra of IRAS 08339+6517 (top) and the companion galaxy 2MASX J08380769+6508579 (bottom). 
The most important emission lines are labeled. The spectra have been scaled down in flux in order to distinguish the faint lines. 
They are not corrected by reddening. Note the high continuum level and the HI absorptions due to the underlying stellar 
population in the spectrum of the main galaxy. 


arm (both showing lowest values). The equivalent width 
of the Ha emission can be used to estimate the age of 
the bursts. We h ave perfo rmed t his es timation using the 
STARBURST99 (iLeitherer et al.l T 999l ) models for an in¬ 
stantaneous starburst (see §4.1 for details). The results for 
each object are found in Tabic 0 

Studying the complex struct ure of the Lyq emission 
observed in IRAS 08339+6517, iMas-Hesse et alJ ll2003|) 
conclude that an expanding superbubble of ionized gas 
powered by the starburst is moving toward the ob¬ 
server at 300 km s -1 with respect to the central Hil 
region. This is supported by the detec tion of X-ray 
em issio n wit h Einstein iGioia et all H990) and ROSAT 
iStevens fe StricklanJl998tT observatories. The strong Ha 
emission detected throughout the disk of the galaxy is 
compatible with this scenario. 

3.4. Intermediate resolution spectra 

The slit position was placed crossing the center of both 
galaxies, with a position angle of 138°. In Figure0we show 
the slit position over the Ha map of IRAS 08339+6517; 
the Ha flux distribution along the slit is shown in Figure^ 
We have also marked in this figure the apertures used 
to extract the individual 1-D spectra. The two brightest 


peaks correspond with the central object (knot #1) and 
knot #2 that we previously detected in the Ha image. We 
have also extracted two external zones, A and B, where 
weak peaks of Ha emission are detected. All four aper¬ 
tures have a size of 1" x 1". But we have also extracted 
four additional spectra with increasing apertures to ana¬ 
lyze the dependence of the properties of the ionized gas 
with the aperture size. These new apertures have a size 
of 2.5", 5", 9", and 12" x 1", and are also indicated in 
Figure |21 The best S/N ratio is found in the largest aper¬ 
ture, corresponding to the integrated spectrum of IRAS 
08339+6517. This aperture collects approximately 96% of 
the total flux of the galaxy observed by our slit. Only 
one spectrum of the companion galaxy was extracted. Its 
aperture size was 6.5" x 1". 

The wavelength- and flux-calibrated spectra of the 
12" x 1" aperture for IRAS 08339+6517 and the single 
aperture for the companion galaxy are shown in Figure 0J 
The spectrum of IRAS 08339+6517 shows a high conti¬ 
nuum level and an intense underlying stellar absorption 
in Hi Balmer lines. All of them, except Ha, show in- 
tense absorption wings . This was previously noted by 
IMas-Hesse et al.l 1 2003 ). Furthermore, we have detected 
the main absorption features present in this spectral 
range: Can H A3933, G-band A4304, Mgi AA5167,87 and 
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Table 4. Dereddened line intensity ratios with respect to I(H/?)=100 for the two galaxies and the selected knots of 
IRAS 08339+6517. Correction for underlying stellar absorption in Hi Balmer lines has been also applied. 


Line 

f(A) 

IR AS 08339+6517 

#1 

#2 

A 

B 

Companion 

3727 [O II] a 

0.27 

288 ± 11 

222 ± 20 

242 ± 17 

322 ± 115 

197 ± 100 

310 ± 44 

3869 [Ne III] 

0.23 

16.8 ± 3.2 

16.7 ± 4.9 

11.5 ± 2.4 



17.6: 

3889 He I + H8 

0.22 

10.6 ± 3.1 

5.3: 

8.5: 




3968 [Ne III] +H7 

0.21 

9.7 ± 0.9 

5.8: 

4.3: 



6.9: 

4101 H<5 

0.18 

25.9 ± 3.1 

25.9 ± 3.6 

25.9 ± 4.6 



23: 

4144 He I 

0.17 

2.4: 






4340 Hy 

0.135 

46.9 ± 3.2 

46.9 ± 3.7 

46.9 ± 3.9 



42: 

4363 [O III] 

0.13 

0.70: 

0.8: 

1.0: 




4471 He I 

0.10 

3.17 ± 0.50 

1.4: 

1.3: 




4658 [Fe III] 

0.05 

1.20 ± 0.55 

0.8: 





4686 He II 

0.05 


1.9: 





4861 H/3 

0.00 

100 ± 3 

100 ± 4 

100 ± 4 

100 ± 30 

100 ± 25 

100 ± 13 

4959 [O III] 

-0.02 

56 + 2 

63 + 4 

49 + 4 

34: 

64 ± 24 

44 + 9 

5007 [O III] 

-0.03 

168 ± 4 

190 ± 10 

163 ± 9 

102 ± 36 

211 ± 62 

135 ± 19 

5200 [N I] 

-0.05 

2.38 ± 0.16 

2.1: 

2.3: 




5516 [Cl III] 

-0.15 

0.6: 






5538 [Cl III] 

-0.16 

0.5: 

0.7: 





5755 [N II] 

-0.21 

0.82 ± 0.29 

1.1: 

0.6: 




5876 He I 

-0.23 

11.6 ± 0.8 

13.3 ± 1.8 

11.7 ± 1.9 



11 + 4 

6300 [O I] 

-0.30 

5.5 ± 0.8 

3.7 ± 1.1 

5.0 ± 1.0 




6548 [N II] 

-0.34 

29.4 ± 1.5 

27.1 ± 2.4 

26.9 ± 2.9 

23 + 9 

26 + 6 

17.6 ± 4.9 

6563 Ha 

-0.34 

286 ± 10 

286 ± 14 

286 ± 15 

286 ± 83 

286 ± 78 

286 ± 37 

6584 [N II] 

-0.34 

79.3 ± 3.1 

78.9 ± 5.6 

76.4 ± 5.5 

62 ± 19 

81 ± 19 

58.2 ± 9.3 

6678 He I 

-0.35 

2.91 ± 0.38 

2.3: 

3.4: 




F(R/3) b 


192 ± 6 

36.5 ± 1.3 

28.7 ± 1.0 

0.329 ± 0.097 

0.324 ± 0.081 

2.39 ± 0.31 

am 


0.22 ± 0.02 

0.30 ± 0.02 

0.29 ± 0.02 

0.18 ± 0.03 

0.25 ± 0.03 

0.18 ± 0.03 

Wabs C 


1.8 ± 0.1 

1.1 ± 0.1 

0.9 ± 0.1 

1.5 ± 0.1 

1.5 ± 0.1 

1.5 ± 0.1 

—W (Ha) c 


104 ± 2 

140 ± 6 

93 + 4 

18 + 2 

26 + 3 

27 ± 0.02 

-W(H/3) C 


19.0 ± 0.8 

25 + 2 

17.6 ± 1.5 

5.6 ± 1.1 

6.7 ± 0.7 

7.5 ± 0.2 

-W(H 7 ) c 


5.9 ± 0.3 

4.9 ± 0.4 

5.7 ± 0.6 




-W(H<5) C 


2.1 ± 0.4 

1.2 ± 0.2 

2.3 ± 0.4 




-W([0 III]) C 


27.2 ± 0.9 

24 ± 1 

29 + 1 

4.0 ± 0.4 

6.7 ± 0.4 

7.9 ± 0.4 


a [O II] AA3726,3729 emission flux. 
b In units of 10 -15 erg s -1 cm -2 . 
c In units of A. 


Na D1 A5889, as well as some Hei absorption. These fea¬ 
tures, together the high continuum level of the spectrum, 
can be interpreted as the product of a substantial pop¬ 
ulation of older stars, with ages much greater than 10 
Myr. Consequently, at least two separate stellar popula¬ 
tions are present in the galaxy (see §4.1). The spectrum of 
the companion object is noisier but shows a rather flat low 
continuum level with no important absorption features, 
characteristic of an H n region-dominated zone. 

Fluxes and equivalent widths of the emission lines were 
measured using the SPLOT routine of the IRAF package. 
This task integrates all the flux in the line between two 
given limits and over a fitted local continuum. All the 
emission lines detect ed were corrected for reddening us¬ 
ing the ISeatonl (Il979il law and the corresponding value of 
C(H/3) derived in the way we explained in §3.1. The H /3 
line flux, F(H[3) and the equivalent width of several lines 


(Ha, H/3, Hy, H<5 and [Oill]) are given in Tabled as well 
as the reddening-corrected line intensity ratios relative to 
H/3 and their uncertainties. Colons indicate uncertainties 
greater than 40%. Note that the values obtained for the 
Ha equivalent width of #1, #2, IRAS 08339+6517 and 
the companion galaxy using the spectra and the Ha im¬ 
age are similar to within the uncertainties. 

3.4.1. Physical conditions of the ionized gas 

We have studied the physical conditions and chemical 
abundances of the ionized gas from our 1-D spectra. 
Although the weak [O ill] A4363 emission line is barely 
detected in the brightest galaxy, the underlying stellar 
absorption affecting 77+ does not permit a proper mea¬ 
surement of that auroral line, except for the 12" aper¬ 
ture (see below). The oxygen abundance was derived then 
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Table 5. Chemical abundances of the galaxies and analyzed knots inside IRAS 08339+6517. 


Object 

IRAS 08339+6517 

#1 

#2 

A 

B 

Companion 

T e [ O III] 3 

8700 

7900 

7900 

8700 

7700 

9050 

T e [0 II] b 

9100 

8500 

8500 

9100 

8400 

9300 

12+log 0/H c 

8.45 ± 0.10 f 

8.55 ± 0.10 

8.53 ± 0.10 

8.42 ± 0.10 

8.58 ± 0.10 

8.38 ± 0.10 

12+log 0/H d 

8.71 ± 0.20 

8.71 ± 0.20 

8.70 ± 0.20 

8.63 ± 0.20 

8.72 ± 0.20 

8.62 ± 0.20 

(0 ++ +0+)/0+ 

1.53 

1.85 

1.63 

1.29 

2.04 

1.81 

log N+/0+ 

-0.94 

-0.89 

-0.94 

-1.10 

-0.86 

-1.13 

12+log N/H 

7.51 

7.65 

7.59 

7.21 

7.72 

7.26 

12+log Ne ++ /H+ 

7.54 

7.73 

7.56 



7.48: 

12+log He+/H+ 

10.91 ± 0.03 

10.98 ± 0.06 

10.92 ± 0.07 



10.90 ± 0.13 

12+log Fe ++ /H+ 

5.95 ± 0.16 

5.91: 





[O/H]/[O/H] 0 e 

0.62 

0.78 

0.74 

0.58 

0.83 

0.52 


a In units of 10 4 K, calculated using empirical calibrations of lPilvuginl (l200lT) _ 

b In units of 10 4 K, T e [OII] calculated from empirical calib ratio n and the relation given hv lGa.rnettl ill 992h 
c Determined using the empirical calibration_of|Pil^iginjj|200lJ) 
d Determined using the [N II]/Ha ratio jDenico^^^aL 12002) 

6 Asuming 12+log[0/H]o = 8.66+0.05 llAsnhin^^^^l2004ll 
f The oxygen abundance derived for this object using the direct method is 8.42 


using the IPilvuginl (l200lll empirical calibration. The elec¬ 
tron temperatures have been estimated from the T([0 ill]) 
and T([Oii]) pairs that reprod uce the tot al oxygen abun¬ 
dance obtained by applying the lPilvurin ; 1200 1 : empirical 
method (see next section). Like Pilyugin, we have assumed 
a two-zone scheme, the usual relation O/H = 0 + /H + + 
0 ++ /H + and the linear relation between T([Oiii]) and 
T([O ill) based on photoionization models obtained by 
lOarnettl l|l992tl . The finally adopted T[Oiii] and T[0n] 
for each burst are shown in Table EJ We have no determi¬ 
nation of electron densities, so we assumed that N e = 100 
cm -3 for all the objects. 

The spectrum obtained using the 12" aperture allows 
a tentative determination of the electron temperature of 
the ionized gas because of the detection of the [O ill] A4363 
and [N 11 ] A5755 emission lines. These weak auroral emis¬ 
sion lines have an error of ~54% and 34%, respectively, 
giving T([0 hi]) = 8900™ and T([Nn]) = 9000™. 
These values are in agreement with those we obtained 
using the empirical calibration. F urthermore, the direct 
determination also agrees with the Garnettl ( 1992 1 linear 
relation between the high and low electron temperatures. 

3.4.2. Abundance Analysis 

IPilvuginl ll200lh performed a detailed analysis of observa¬ 
tional data combined with photoionization models to ob¬ 
tain empirical calibrations for the oxygen abundance from 
the relative intensities of strong optical lines. He used both 
the i ?23 and P (an indicator of the hardness of the ionizing 
radiation) parameters for his calibration, the so-called P- 
method , which can be used in moderately high-metallicity 
Hu regions (12 +log O/H > 8.15). We have used this em¬ 
pirical calibration to derive the total oxygen abundance 
in IRAS 08339+6517 and its companion galaxy. We have 


also esti mated the to tal oxygen abundance making use 
of the iDenicold et, nil [ 2002 1 empirical calibration, which 
involves the [N 11 ] A6583 /Ha ratio. The oxygen abun¬ 
dances obtained using both methods are shown in Tabled 
However, the O/H ratios obtained using the last calibra¬ 
tion seem to be sys tematically higher than those obtained 
using thelPilyuginl l200ll) method, as we r eported in a pre¬ 
vious study ijLopez-Sanchez et al ] 2004bl) . Here, an offset 
of 0.15-0.20 dex is found betw een both calibrations. We 
prefer to consider the IPilvuginl ll200lll abundances as the 
more appropriate ones because his calibration take into 
account the excitation parameter of the ionizing radiation. 

All the ionic abundances shown in Table [5] (except 
He + and Fe ++ ) have been calculated using the IRAF 
NEBULAR task from the intensity of collisionally ex¬ 
cited lines, assuming the T([Om]) for the high ioniza¬ 
tion potential ions 0 ++ and Ne ++ , and T([Oii]) for the 
low ionization potential ions 0 + and N + . We have as- 
su med the st a ndard ioniza tion correction factor (icf) by 
IPeimbert fc Costerol (. 1969 1. N/O = N+/0+, to derive the 
nitrogen abundance of the objects. The He + /H + ratio has 
been derived from the Hei A5876 line observed in each 
knot, usin g the p redicted line emissivities calculated by 
ISmith et alJ dl996lh We have also correct ed for collisional 
contri bution following the calculations bv lBeniamin et ’all 
d2002lh Self-absorption effects were not considered for the 
determination of He + /H + but only errors derived from 
the observed uncertainties in line measurement. 


We have detected the [Fein] A4658 emission line in 
the integrated spectrum of the main galaxy and in knot 
#1. The calculations for Fe ++ have been done with a 34 
l evel model-atom that uses the collision strengths of 
( 1996 1 and the transition probabilitie s of Quineti 


Applying the relation given bv lRodriguez fc Rubinl 


Zhang 

1996. 

2005). 
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Fig. 5. Position-velocity diagram for the slit position ob¬ 
served, analyzed in 1.1" bins. The horizontal bars represent 
the uncertainty of the Gaussian fitting for each point. Note 
that the vertical axis has been broken for clarity. Northwest is 
up. 


Fig. 6. Position-velocity diagram shown in FigureElcom pare d 
with that obtained using the Hi data from iGannon et al. I 
(1200415 (open circles). The outflow of the neutral gas is clearly 
observed, as well as the decoupling of the ionized gas at the 
NW arm of IRAS 08339+6517. 


we derived a total Fe abundance of 12+log(Fe/H)=6.11 
and log(Fe/0)=—2.60 for the galaxy. 

In general, we find that the O/H ratio of the knots 
inside IRAS 08339+6517 are similar [around 8.5, in units 
of 12 + log(0/H)], except object A, which has a slightly 
lower value [12+log(0/H)=8.42]. On the other hand, the 
oxygen abundance of the companion galaxy is a bit lower 
than that derived for IRAS 08339+6517. The N + /0 + ra¬ 
tios are similar for all the objects, within the expected 
values for spiral and irregular galaxi es with the same oxy¬ 
gen abundance dOarnett et al.lll99^1 . These facts suggest 
that both galaxies have suffered a comparable degree of 
chemical enrichment despite their different absolute mag¬ 
nitudes. 

3.4.3. Kinematics of IRAS 08339+6517 

The kinematics of the ionized gas were studied via the spa¬ 
tially resolved analysis of the Ha line profile along the slit 
position. We have extracted zones 6 pixels wide (1.1" long) 
covering all the extension of the line-emissio n zones . The 
Starlink DIPSO software llHowarth fc Murravlll99 Q.) was 
used to perform a Gaussian fit to the Ha profiles. For each 
single or multiple Gaussian fit, DIPSO gives the fit param¬ 
eters (radial velocity centroid, Gaussian sigma, FWHM, 
etc) and their associated statistical errors. In Figure0we 
show the position-velocity diagram obtained, indicating 
the position of the studied objects. All the velocities are 
referred to the heliocentric velocity of the center of IRAS 
08339+6517 (5775 km s” 1 ). 

The main galaxy has an apparent global rotation pat¬ 
tern, but it also shows local velocity fluctuations that 
could be produced by the effects of the very recent star 
formation activity (gas motions due to the combination 
of the winds of massive stars and supernova explosions) 


or by distortions in the gas associated with interaction 
effects (with the companion galaxy or, perhaps, with an 
hypothetical pre-existing external object that could have 
merged with IRAS 08339+6517 in the past). Interaction 
features are clearly detected in the outer northwest area of 
IRAS 08339+6517, where the ionized gas does not follow 
the general kinematics but presents a continuous variation 
between +45 and +130 km s -1 . These positions coincide 
with the Ha emission detected in knot B and the NW 
arm. The companion galaxy seems to shows an apparent 
rotation behavior in its southeast area (its brighter zone), 
although an important velocity gradient (around +200 km 
s” 1 ) is found in the side opposites the main galaxy. This 
fact suggests again that gravitational interactions have af¬ 
fected the dynamics of the ionized gas. 

Assuming the global kinematics of the galaxies as cir¬ 
cular rotation, we can make a rough estimate of their 
Keplerian mass, Mk ep , taking into account the half of 
the maximum velocity difference in the galaxy and the 
half of the spatial separation corresponding to these max¬ 
ima. We do not consider in this calculus the effect of the 
velocity gradients presented in both galaxies. We obtain 
masses of (10±3)xl0 9 M© and (8±2)xlO s M© for IRAS 
08339+6517 and its companion galaxy, respectively, as¬ 
suming circular orbits, Keplerian dynamics and an incli¬ 
nation of 90°. We remark again that these values are only 
very tentative estimations of the masses of the galaxies. 

We have compared the Keplerian mass, the neutral 
gas mass, Mrj, and the ionized hydrogen mass, M^ii 
(see Table 0 of the both galaxies and find substantial 
differences between them. IRAS 08339+6517 has higher 
Mhii/Mhi and Mk ep /M r i ratios (0.013 and 9.1, re¬ 
spectively), while the companion object has 0.0004 and 
1.1. However, if we co nsider that the tidal Hi gas found 
between both objects ijCannon et al. Il2004 ). with a mass 
of Mhi = (3.8±0.5)xl0 9 M©, have been expelled from 
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the main galaxy, the new Mk ep /M r i ratio of IRAS 
08339+6517 is 2.0, similar to that found in the companion 
galaxy. 


The light-to-mass ratios, L B / Mk ep and Lb/M m, are 
also very different for both galaxies. The main object has 
Lb/M hi ~ 61, much higher than the typical values for 
spiral and ir regula r galaxies, that are between 1.6 and 
11.2 dBettoni et, a,l.lEoQ3 1. The companion galaxy has a 
value between those limits (Lb/Mhi ~ 4.2), similar to 
an Sbc galaxy. If we again consider that the tidal neutral 
gas belongs to IRAS 08339+6517, the new light to neu¬ 
tral gas ratio is ~13.7, only a bit higher than the expected 
value for an Sa galaxy (11.2 ± 2.0). The light-to-cold dust 
ratio, L B /M dust ~ 1.5 x 10 4 , is also the expected one 
for young spirals dBettoni et ~afl l2003ih The mass of the 


dust was calcula ted using the FIR 60/im and 100/im flux 
densities jjMgfjhir^ a! U 1991^ tog ether with the relation 


given bv f^uch^eier et a,1.1 



(1 


(their last equation in 


p. 678). All these facts support the idea that the tidal Hi 
gas has been formed from material mainly stripped from 


IRAS 08339+6517. 


In order to improve our knowledge of the kinematics of 
the system, we have compared the position-velocity dia¬ 
gram obtained from our analysis of the Ha line profile in 
the optical spectrum with the position-velocity diagram 
ob tained from HI data. We have used the Figure 2 of 
[Cannon et~al~l J2004 to extract the velocity of the neu¬ 
tral gas along our slit position (with P.A. 138°) in 10" 
steps. We referred all the Hi velocities to the heliocen¬ 
tric velocity of the center of IRAS 08339+6517 (5775 km 
s _1 ). In Figure El we show the position-velocity diagram 
combining HI and Ha data. 

The general correspondence between both diagrams 
is rather good: the centers of the two galaxies and the 
kinematics of IRAS 08339+6517 p ractically coinc ide in 
both neutral and ionized gas. ICaxmon et all ( 2004 1 have 
remarked that, although the companion galaxy shows a 
component of solid-body rotation in neutral gas, the signs 
of rotation are less prominent in the main galaxy. But it 
could be masked in part by the outflow of the neutral gas 
in the direction of the observer, as is clearly observed in 
Fi gure El This ou tfl ow of ne utral gas was firstly suggested 
bv IGonzalez-Delgado et aT ■ •199^1 because of the excess 
absorption detected in the blue wing of the Ovi A1032 
profile. They also noted that the broadening of the Lya, 
Sill A1260 and Cll A1335 absorption lines (the metal¬ 
lic lines are also resolved in two components) indicated 
large-scale motions of the interstella r gas, with a velocity 
even higher than 3000 km s -1 . IGonzalez-Delgado et all 
: 1998 1 suggested that the kinematics of the warm ion¬ 
ized interstellar medium are not driven by their gravi¬ 
tational potential but by the dynamical consequences of 
the violent star form ation processes taking place in the 
galaxy. iKunth et~ahl dl998l) attributed the existence of a 
secondary Lya emission peak to the chaotic structure of 
the interstellar medium. Our Ha data suggest that, al¬ 
though these significant distorsions are found in the ion¬ 


ized gas, the kinematics of the galaxy still have an appar¬ 
ent important component of solid-body rotation. 

However, we find a clear divergence between the kine¬ 
matics of the neutral and the ionized gas at the NW arm 
of IRAS 08339+6517, where object B is located: while the 
first has a velocity of —50 km s -1 with respect to the 
center of the galaxy, the second shows a velocity around 
+100 km s -1 . Thus, the ionized gas seems to be decoupled 
from the expanding motion of the tidal tail of neutral gas. 
This suggests that the interaction that has created the HI 
tail is more complex that we thought. Knot B is specially 
interesting because it has the hightest oxygen abundance 
(12 + log O/H = 8.58), being the higher of all the spec¬ 
troscopically observed regions (see Table 0 . 

4. Discussion 

4.1. Ages of the bursts and stellar populations 

4.1.1. Ages from VF(Ha) 

One of the best methods of determining the age of the 
recent bursts of star formation is through the Ha equiva¬ 
lent width, W(Ha), since i t decreases with time. We have 
used the STARBURST99 libeitherer et a,hill 9991 spectral 
synthesis models to estimate the age of each knot, compar¬ 
ing the predicted Ha equivalent widths with the observed 
values. Spectral synthesis models with two different metal- 
licities were chosen (Z/Z 0 = 1 and 0.4, the appropriate 
range of metallicities for the objects, see Table 0 , both 
assuming an instantaneous burst with a Salpeter IMF, a 
total mass of 10 6 M 0 and a 100 M 0 upper stellar mass. 
The ages obtained are compiled in Table 0 All of them 
show a very young age for the last star formation episode, 
between 3.5 and 6 Myr, being the center of the galaxy 
(knot #1) the youngest one. We can adopt an average age 
of 4 - 5 Myr for the most recent star formation episode in 
IRAS 08339+6517 and 5-6 Myr in its companion galaxy. 

4.1.2. Ages from evolutionary synthesis models 

iGonzalez-Delgado et al . (Il998h used evolutionary synthe¬ 
sis models to constrain the star formation history of IRAS 
08339+6517 from their UV spectra. They found that a 6- 
7 Myr burst with M up > 30 M 0 or 9 Myr continuous 
star formation with M up < 30 — 40 M 0 are compatible 
with their observations. They noted that the weakness of 
the stellar UV absorption lines and the small Ha equiv¬ 
alent width suggested significant dilution by an underly¬ 
ing stellar population. Our estimated ages from W{ Ha) 
are a little lower than those they predicted. These au¬ 
thors also compared the C IV and Si IV profiles with those 
of the starburst galaxy NGC 1741, suggesting that IRAS 
08339+6517 is probably in a more advanced evolutionary 
stage. But the weakness of these lines could also be ex¬ 
plained if the IMF is steeper or there exists a dilution 
of the stellar lines due to the underlying population. The 
latest study of NGC 1741 (memb er C in HCG 31) per¬ 
formed bv iLopez-Sanchez et al . ll2004ah estimates an age 
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Fig. 7. Co lor-c olor diagrams showing the STARBURST99 
Lejthereii_eGalJll£99l) (blu e/con tinuous line) and PEGASE.2 
Fioc ,+ Rocca.-Volmera.ngFll 997f) (dashed/red line) models for 
an instantaneous burst with a Salpeter IMF and Z/Zq = 0.4 
compared with our observed values for IRAS 08339+6517 (cir¬ 
cle) and its companion galaxy (triangle). We have also in¬ 
cluded the derived values for the underlying component in 
IRAS 08339+6517 (cross) and the burst component (square). 
We indicate some age ticks (in Myr) along the plotted evolu¬ 
tionary tracks for the STARBURST99 model. 



Fig. 8 . F [Q hi] A5007 versus W (H/3) using the models by 
IStasinska et al .' 11200 1+ Tracks correspond to sequences of dif¬ 
ferent metallicities and electronic densities. Each symbol marks 
the position of the models at 1 Myr interval, starting in the 
upper-right corner of the diagram with an age of 1 Myr. 


of around 5 Myr and E(B—V) =0.06 for this galaxy. Thus, 
IRAS 08339+6517 and NGC 1741 have very similar ages, 
but the first has a larger reddening contribution. Perhaps 
the effects of an underlying population is the most impor¬ 
tant factor in the dilu ti on of the UV stellar lines in the 
iGonzalez-Delgado et al.1 ( 1998 ) study. 

We have comp ared our reddening-cor rected colors with 
STAR BURST9 9 dbeitherer et al l il99Q') and PEGASE.2 
( Fioc fc R.occa-Volmerangel 1199711 models to get an ad¬ 
ditional estimation of the age of the two galaxies. We 
have chosen these two models because while the firsts are 
based on Geneva tra cks, the seconds use Padua isochrones 
iBertilli et al.1 Il99~4 ) in which thermally pulsing asymp¬ 
totic giant branch (TP-AGB) phases are included. We 
assumed an instantaneous burst with a Salpeter IMF, a 
total mass of 10 6 Mq and a metallicity of ZIZq = 0.4 
and 1 for both models. In Figure [7| we plot four different 
color-color diagrams comparing both STARBURST99 and 
PEGASE.2 ZIZq = 0.4 models with the observed values. 
We find a good correspondence using our optical values al¬ 
though PEGASE.2 give older ages than STARBURST99. 
In particular, the V — R color of IRAS 08339+6517 im¬ 
plies an age of 40 - 60 Myr using STARBURST99 models, 
but 450 - 600 Myr when it is compared with PEGASE.2 
models. Furthermore, the NIR colors derived from 2MASS 
(reddening-corrected, see Table 13 are not in agreement 
with optical colors, specially in the case of the companion 
galaxy. The optical colors suggest an age of ~ 30 - 50 Myr 
for IRAS 08339+6517 and older than 200 Myr for its com¬ 
panion. The NIR colors imply ages older than 1-2 Gyr 
for both systems. The considerable difference in the age 
derived for the companion galaxy using the W (Ha) and 


the optical and NIR colors indicates that the old popula¬ 
tions practically dominate the observed flux at these long 
wavelengths. 

The fact that the ages derived from colors are older 
than those obtained from W(Ra) confirms that the ef¬ 
fects of the underlying old population are important. The 
contributi on of the old ste l lar p opulation can be checked 
using the IStasinska et al.l ( 200 1) models of HII regions 
ionized by an evolving starburst embedded in a gas cloud 
of the same metallicity. We have chosen two models, both 
with metallicity Z/Z 0 = 0.25 and 1 and with a total mass 
of 10 3 and 10 6 Mq. In Figure 0 we plot our observational 
values of W (H/3) and the [O ill] A5007/H/3 emission line 
flux (see Table and compare them with the theoreti¬ 
cal models. We find that all are located in the left area 
of the diagram, indicating a strong contribution of the 
underlying stellar continuum. Objects #1 and show 
the best correspondence because they are regions mainly 
dominated by the recent star formation activity. The age 
derived for the objects following Figure[B]is between 5 and 
6 Myr, in agreement with those derived from IT(Ha). 

4.1.3. Analysis of the surface brightness profiles 

In order to investigate the importance of older stellar 
populations we have performed an analysis of the surface 
brightness profiles of the two galaxies. We have taken con¬ 
centric surfaces at different radii from the center of each 
system and calculated the integrated flux inside each cir¬ 
cle of area A (in units of arcsec 2 ). Then, the mean sur¬ 
face brightness inside this circle, SBx (in units of mag 
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Table 6. Structural parameters of IRAS 08339+6517 and its companion galaxy. The total luminosity in each band 
and the relative contribution of the underlying component at some radii are also shown. 


Filter 

Mo 

(mag ”- 2 ) 

a 

(kpc) 

IRAS 08339+6517 
Ltotal 

(10 43 erg s -1 ) 

% UC 
r = 2.5” 

% UC 
6” 

% UC 
14” 

Mo 

(mag ”“ 2 ) 

a 

(kpc) 

Companion 

Ltotal 

(10 41 erg s -1 ) 

% UC 
r = 3.5" 

t/ a 

21.73 

4.49 

2.15+ 0.08 

1 

13 

92 

19.41 

0.76 

7.1 ± 0.5 

84 

B 

21.77 

4.31 

3.09 ± 0.11 

2 

17 

99 

19.51 

0.73 

14.0 ± 0.6 

92 

V 

21.61 

4.31 

1.72 ± 0.07 

3 

16 

96 

19.31 

0.74 

9.2 ± 0.5 

95 

R 

21.28 

4.21 

2.20 ± 0.09 b 

3 

14 

96 

19.26 

0.83 

11.9 ± 0.6 

95 


a The errors in the fitting parameters in U filter are higher than those found in the other filters. See text. 
b Considering the Ho emission (see Table 0, the total luminosity in R filter is (2.08 ± 0.10) x 10 43 erg s' 



Fig. 9. Surface brightness and color profiles for IRAS 
08339+6517. The dashed lines in the surface brightness dia¬ 
gram represent exponential law fits to the profiles. The red 
horizontal lines in the color profile diagrams indicate the av¬ 
erage color derived for each system. The dashed vertical lines 
indicate the radius of the seeing (0.4”) and the interval used 
to perform the fit (between 14” and 22”). The peak at 24” is 
consequence of the bright star located at the SW of the galaxy. 

arcsec -2 ), is derived using the relation: 


Fig. 10. Surface brightness and color profiles for the compan¬ 
ion galaxy. The dashed lines in the surface brightness diagram 
represent exponential law fits to the profiles. The red horizon¬ 
tal lines in the color profile diagrams indicate the average color 
derived for each system. The dashed vertical lines indicate the 
radius of the seeing (0.4”), and the interval used to perform 
the fit (between 2” and 7”, the apparent optical radius of the 
galaxy derived from the deep optical images). Note the bump 
at ~3.5” in 17-filter, it is indicating the star-formation zone 
off-center this dwarf galaxy. 


SB X = mx + 2.5 log A, 


( 2 ) 


being mx the magnitude in the filter A'. The surface 
brightness, mx, is the flux per square arcsec in the ring 
defined by two successive apertures. This simple technique 
is not adequate for the study of o b jects pr esenti ng irregu¬ 
lar or complex morphologies ( Cairos et aflbOQll) . but it is 
valid for approximately circular compact objects like IRAS 
08339+6517 and its companion galaxy. In Figure 0 we 
show the surface brightness for U, B , V and R filters, mj, 
Hb, Mv and Hr, versus its radius, for IRAS 08339+6517, 
whereas the surface brightness profiles for the companion 
galaxy are shown in Figure m In both Figures 151 and [TUI 
we also include the radial color profiles (U — B), {B — V) 
and (V — R) derived by direct subtraction of light profiles. 

In the main galaxy, the surface brightness profiles can 
be separated into two structures, indicating the existence 


of a low surface brightness component underlying the star- 
burst that could be attributed to a disk-like structure in¬ 
side the galaxy. We have performed an exponential law 
fitting to the profiles, following the expression: 

I = I 0 exp(-r/a), (3) 

which describes a typical disk structure: I 0 is the central 
intensity and a is the scale length. The fits were performed 
between r m i n = 14” and r max = 22” and are plotted over 
each profile with a dashed line. The upper limit was chosen 
in order to avoid contamination by the bright star at the 
SW of the galaxy (see Figure [IJ; its effect on the surface 
brightness can be observed as the peak at a radius of 24". 
Thus, it introduces an additional uncertainty to the linear 
fits to the underlying stellar emission. The fitting struc¬ 
tural parameters are indicated in Table Remark that an 
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exponential intensity profile is typical not only for disks 
but also for low-mass spheroids, so it should not be taken 
as evidence that the 3D geometry of IRAS 08339+6517 is 
approximated best by a disk. Note that the fit for the U fil¬ 
ter has highest error because the image is not deep enough 
to reach the faint surface brightness level needed to per¬ 
form a good analysis. The fits can be used to determine 
the relative contribution of the underlying component and 
the burst (we will designate as burst the luminosity that 
can not be explained by the underlying component) to the 
luminosity at different radii. The relative contribution of 
the underlying component at several radii is also indicated 
in Tabic |6] it is only 2% at the central areas but equal to 
the contribution of the burst at ~10". From 11" it domi¬ 
nates the total luminosity of the galaxy. The burst to total 
luminosity ratio is ~85%. 

The variations of the radial color profiles also indi¬ 
cate the different stellar populations present in the main 
galaxy. Its central areas show bluer colors, but for radii 
> 14" the color profiles have no gradients and show va/- 
lues redder to the ones derived from the integrated pho¬ 
tometry. Surprising, we note a central color excess in the 
U — B profile inside radius ~ 2". It does not seem to be ar¬ 
tificial because its size is larger than the seeing (~ 0.8"). 
In Figure we present an U — R color map of IRAS 
08339+6517, as well as the U — R color profile. The color 
map shows a ring-structure with bluer colors, that includes 
the two central bright bursts (see Figure El, around the 
nucleus of the galaxy. The U — R color map also repre¬ 
sents the dust distribution throughout the galaxy, that is 
somewhat inhomogeneous, as we previously note in §3.4.3. 
The nucleus of the galaxy seems to be dustier than its sur¬ 
roundings, being probably the explanation of the central 
color excess observed in both U — B (Figure |9j) and U — R 
(Figure El color profiles. This feature is not observed in 
the others color profiles because extinction in U filter is 
much more effective. 

We have estimated the colors of the underlying com¬ 
ponent and the burst using the fits to the light profiles 
derived above. They are plotted in Figure 0 As we ex¬ 
pected, the underlying component has significantly redder 
colors than the burst, supporting that it is composed by 
an older population. Their U — B ~ —0.08 and B — V = 
0.16 colors imply a minimal age between 200 and 300 Myr, 
although it rises to 700 - 900 Myr considering the V — R 
= 0.30 color. For the burst component we find U — B = 
—0.51, B — V = 0.02 and V — R = 0.23, suggesting an age 
between 15 and 25 Myr using both STARBURST99 and 
PEGASE.2 models. These values are not as young as those 
obtained from W(Hq!); it may indicate that, besides the 
disk-like underlying structure, which has been removed, 
an additional contribution of older stars also exists in the 
central areas of the galaxy. 

On the other hand, the fits in the companion galaxy 
were performed between r m i n = 2" and r ma x = 7". They 
do not show two structural components, although devi¬ 
ations from the fits can be observed in the inner area 
and at r ~ 3.5" (see Figure m where the starburst is 



-0.25 0.00 0.25 0.50 



Fig. 11. U — R color map of IRAS 08339+6517, revealing a 
probable inhomogeneous distribution of the dust in the inner 
areas of the galaxy. Remark the blue ring-like structure around 
the nucleus of the galaxy. We also show the U — R color profile 
for IRAS 08339+6517. Note the color excess for radii < 1". The 
red horizontal line indicate the average color, U — R — —0.25. 


located and accounts for ~ 10% of the surface bright¬ 
ness. It seems to be somewhat more evident in the pro¬ 
file derived for the C/-filter. In general, the contribution 
of the starburst is diluted under the old stellar compo¬ 
nent, that accounts of 99% of the total luminosity of this 
d warf galaxy. High spatial resolution photometry in NIR 
dNoeske et all2003h should be performed for a proper sep¬ 
aration of these two stellar components. The radial color 
profiles are rather constants, with values similar to those 
derived from aperture photometry, also supporting the 
idea that the galaxy is mainly dominated by this old stel¬ 
lar population. Remember that the U filter is not deep 
enough to reach a faint surface brightness level; its uncer¬ 
tainty is high for radii greatest than r ~ 6". 


4.1.4. Ages from hydrogen stellar absorption lines 


The age of the stellar populations in a starburst 
galaxy can be estimated analyzing the equiva¬ 
lent widths of hydro gen and helium absorption 


lines. Gonzajez-Delgad o&Mjeithererl ( ll99 9?) and 


iGonzalez-Delgado et al ■ ill 999(1 


presented a synthetic 
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Fig. 12. (Right) H 7 and H<5 emission lines of our 12" spectrum 
(black line) and a fit (long dashed line) using a narrow com¬ 
ponent ontop a broader absorption component (dash-dotted 
line). (Left) Variation of the equivalent widths of the H 7 and 
H<5 absorption lines with the age of the burst, following the 
evo l utionary stel lar population synthesis models presented by 
iGonzalez-Delgado et al.1 ( 1999 1. The dashed horizontal lines 
represent the equivalent widths derived from our spectrum. 
Note the decrement of the equivalent width for t >500 Myr in 
the models. 


grid of stellar Hi Balmer and Hei absorption lines in 
this kind of galaxies and developed evolutionary stellar 
population synthesis models for instantaneous bursts 
with ages between 1 Myr and 1 Gyr, assuming a Salpeter 
IMF between 1 and 80 M© and solar metallicity. The 
models indicate an increment in the equivalent widths 
with age. This analysis has been useful in similar studies 
of y oung star bursts or blue compact dwarf galaxies (e.g. 
iGuseva et al • 2001 ). 

W e have used DIPSO software 1 Howarth fc Murray! 
1990 ) to perform a Gaussian fitting to the H 7 and H4 pro¬ 
files of our 12" spectrum using a narrow Gaussian emis¬ 
sion line profile ontop a broader absorption component 
to estimate their equivalent absorption widths. The fits 
are shown in Figure El We have obtained W a b s ( Hy)= 
3.7 ± 0.1 A and H 4h.riH(5)= 3.9 ± 0. 1 A Comp aring 
these values with the lGonzalez-Delgado et al • ‘TP?)!): mod¬ 
els, we estimate an age between 4 and 7 Myr for the 
stellar population. This age is surprising similar to the 
one derived from the Ha emission and not to that found 
from the analysis of the surface brightness profiles (> 100 
Myr). If the case of continuous star formation is con¬ 
sidered instead of an instantaneous burst scenario, the 
values derived for W a b s ( H 7 ) and W a b s (H5) will imply a 
large r age, ~15 Myr, but sti l l you nger that 100 - 200 
Myr. IGonzalez-Delgado et al.l (1999) models also predict 
than the equivalent widths of Balmer lines in absorption 
decrease for ~500 Myr, so the derived values might also 
imply ages >1 Gyr. It would be in agreement with the es¬ 
timation given by NIR colors. The analysis of H 7 and H(5 
equivalent widths is not useful, in this case, to constraint 
the age of the underlying stellar population because the 


result merely is a luminosity-weighted age for the mixture 
of stars in the central part of the galaxy. At the center 
of IRAS 08339+6517, the light is dominated by the young 
stars recently formed in the starburst, as we see both from 
the Ha image and the fact that it accounts for ~85% of 
the total luminosity at this radius (the contribution of the 
underlying component is r^i 15% at r = 6 ", see Table| 6 l). In 
this sense, an appropriate study of the underlying stellar 
population using absorption lines would require to mea¬ 
sure them at radii r> 15", but our spectra are not able to 
reach these values. 


4.1.5. Ages from the spectral energy distribution 


Finally, we have employed spectral energy distributions 
(SED) to constraint the ages of the stellar populations. 
Although this method is dependent of the interstellar ex¬ 
tinction, our estimation of the reddening contribution us¬ 
ing the Balmer decrement (see §3.1) lets us to perform it 
avoiding the degeneracy problem between reddening and 
the ages of stella r po pul ations. We have made us e of the 
PEGASE.2 code llFioc fc Rocca-Volmerangell997t) to pro¬ 
duce a grid of theoretical SEDs for an instantaneous burst 
of star formation and ages between 0 and 10 Gyr, assum¬ 
ing a Zq metallicity and a Salpeter IMF with lower and 
upper mass limits of 0.1 M 0 and 120 M 0 . Although the 
grid include the ionized gas emission, we have neglected it 
because its contribution to the continuum is rather weak. 
In Figure ED we show our 12 " extinction-corrected spec¬ 
trum and the synthetic continuum spectral energy dis¬ 
tributions derived assuming young (6 Myr) and old (140 
Myr) ages. None of the individual synthetic spectra fitted 
our observed SED, so we constructed a model than com¬ 
bines 85% of the 6 Myr model and 15% of the 140 Myr 
model. We chosen these values following the relative con¬ 
tribution of the burst and the underlying component that 
we derived in the analysis of the surface brightness profiles 
at 6 " (see Table EJ - We note that this combined model is 
in excellent agreement with the shape of our derreddened 
spectrum. However, we want to remark the degeneration 
existing in this kind of studies: several combinations of dif¬ 
ferent relative contributions for old and young population 
models might also explain the observed spectrum. For ex¬ 
ample, a model which combines a contribution of 50% for 
both 6 Myr and 100 Myr models also agrees with it. 


The agreement between observed and synthetic spec¬ 
trum is also found using the Hi Bal mer an d Heiab- 
sorption lines derived from lGonzalez-Delgado et all ( 1999 ) 
models. In Figure 1141 we plot our derredened normalized 
spectrum in the 3700 - 4450 A range compared with their 
Zq models with ages of 4 and 200 Myr. The spectral reso¬ 
lution of the models were degraded to that of our observed 
spectrum. Again, a good fit is found when a combined 
model with a contribution of 15% for the old population 
and 85% for the young population is considered. 


In conclusion, the young stellar population in IRAS 
08339+6516 (the most recent starburst) has an age be- 
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Fig. 13. Spectrum of IRAS 08339+6517 compared with syn¬ 
thetic continuum spe ctral energ y distributions obta ined us¬ 
ing the PEGASE.2 ilFioc fc Ti occa,-Volmera.n gel il997i l code. 
The dotted line represents the observed spectrum uncorrected 
for extinction, whereas the gray/green continuous line is the 
extinction-corrected spectrum assuming C(H/3)=0.22. The up¬ 
per continuous line corresponds to a model with an age of 6 
Myr (young population model), whereas the lower continuous 
one is a 140 Myr model (old population model). The shape 
of our observed derredened spectrum fits with a model with 
a contribution of 85% for the young population and 15% for 
the old population is considered (continuous black line over the 
galaxy spectrum). 


tween 4 and 6 Myr, and it is located in the inner areas of 
the galaxy (radii < 6"). It is superposed to a more evolved 
stellar population with an age not younger than 100 - 200 
Myr, that fits a disk-like profile, and probably formed in 
early bursts. The NIR colors suggest that an older stellar 
population, with ages larger than 1-2 Gyr, also exists 
in IRAS 08339+6517. Its companion galaxy is practically 
dominated by an old population with age > 250 Myr, al¬ 
though a recent starburst of around 6 Myr is found in its 
external areas. 

4.2. The star formation rate 

We have derived the star formation rate (SFR) of IRAS 
08229+6517 using empirical calibrations in different wave¬ 
bands. The SFR is a key parameter for characterizing the 
formation and evolution of the galaxies. However, even in 
the well-observed galaxies of the local Universe, it remains 
quite uncertain because different methods yield different 
values. Much of this uncertainty is related to the unknown 
contribution of the dust obscuration in and around the 
star-forming regions. 

One of the most extended methods of deriving the SFR 
is the use of hydrogen recombination line fluxes, especially 
the Ha flux. Since the flux in a hydrogen recombination 
line is proportional to the number of ionizing photons pro¬ 
duced by the hot stars (which is also proportional to their 
birthrate), the SFR can be easily derived. Th e most re- 
cent calibration for starbursts is that derived hv lKennicnttl 


Fig. 14. Normalized dereddened spectrum of IRAS 
08339+6517 (gray/gree n con tinuous line) compared with 
iGon zalez-Pel gado et aTJ ( 1999 1 models with 4 Myr (dotted 
line) and 200 Myr (dashed line) at Zq metallicity. The best 
fit corresponds to a model with a contribution of 85% for the 
young population and 15% for the old population (continuous 
black line). Note that the models only give values for the 
position of the absorption lines and not for all wavelengths, 
being the reason of the straight lines presented between 3930 
and 4000 A and 4150 and 4300 A. The straight lines are 
artifacts connecting gaps in the dataset of the models. 


( 1998 ). We have used this calibration to estimate the SFR 
in IRAS 08339+6517 from our Ha image, as we previously 
indicated in §3.3, obtaining SFR# a = 9.5 Mq yr -1 (see 
Table 0. Note that we have corrected the Ha flux for both 
reddening and [N n] contamination. 

Many of the problems found in deriving the SFR 
from optical data can be avoided by measuring the far- 
infrared (FIR) and sub-millimeter spectral energy distri¬ 
bution. These are determined by the reradiation by the 
dust of energy absorbed in the visible and UV regions 
of the spectrum. Assuming that the dust completely sur¬ 
rounds the star forming regions, it acts as a bolometer re¬ 
processing the lu minosity pr oduced by the stars. We have 
applied the lKennicuttl (1998) correlation between the SFR 
and the FIR flux using the IRAS satellite data f or IRAS 
08339+6517, /6o=5.90 Jy and /ioo=6.50 Jy llMoshir et alJ 
[m3), deriving SFR^j# = 9.7 Mq yr 1 . This value is in 
excellent agreement with the estimation derived from the 
Ha luminosity, indicating that the value for the extinc¬ 
tion we have adopted seems to be fairly appropriate. We 
have also used others correl ations between SFR and 60/im 
( Condonlll992tl and 15^m jR.oussel et luminosi¬ 
ties (the last one assuming Li^^ m ~ £ 12 /im), finding 

similar values (see Table 0. 

The radio continuum flux can be also used as a star for¬ 
mation indicator. Nearly all of the radio luminosity from 
galaxies without a significant active galaxy nuclei (AGN) 
can b e traced to recent ly formed massive (M > 8 M@) 
stars llGondon et al.lll992tl . Ten percent of the continuum 
emission at 1.4 GHz is due to free-free emission from ex¬ 
tremely massive main-sequence stars (thermal emission) 
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Table 7. Star formation rates (SFR) derived for IRAS 08339+6517 using different relations (see §4.2). All the lumi¬ 
nosities were calculated assuming a distance of 80 Mpc for the galaxy. 


Range 

Luminosity 

Units 

Value 

SFR (M 0 yr" 1 ) 

Calibration 

X-ray 

Lx (0.2 —2.0 keV) 

erg s' 1 

2.81 x 10 41 

61.8 

Ranalli et al. (2003) 

Optical 

L Ho, 

erg s' 1 

(12.0 ± 0 .6) x 10 41 

9.5 ± 0.5 

Kennicutt Q92&) 


L B 

L© 

(6.61 ± 0.24) x 10 10 

1.92 ± 0.07 

Gallagher et al. (1984) 


Lb,uc 

L 0 

(8.95 ± 0.30) x 10 9 

1.1 ± 0.1 

Calzetti (2001) 


L[on] 

erg s 1 

1.21 x 10 42 

13 

K^nnicj^ '(19981 





8.9 

Kewley et al. (2004) 

Far-Infrared 

L FIR 

L© 

5.63 x 10 10 

9.7 

Kennicutt (1998J 


Ll2/xm 

L© 

1.78 x 10 10 

11.6 

Roussel et al. (2001) 


L60 |im 

L© 

3.90 x 10 10 

9.2 

^ondon (l()92l 

Radio 

L1.4 GHz 

W Hz" 1 

2.57 x 10 22 

6.4 

Condon et al. (20021 


and almost 90% is synchrotron radiation from relativistic 
electrons accelerated in the remmants of core-collapse su¬ 
pernovae (non-thermal emission). As the stars that con¬ 
tribute significantly to the radio emission have lifetimes 
r < 3 x 10 7 yr and the relativistic electrons have life¬ 
times r < 10 8 yr, the current radio luminosity is nearly 
proportional to the r ate of massive star f ormation during 
the past r < 10 8 yr I Condon et ~Zl l 2002 ll . 

We have estimated the SFR fro m the 1.4 GHz lumino- 
sity for IRAS 08339+6517 using the lCondon et alJ d2002l) 
correlation. The 1.4 GHz l uminosity was deriv ed from 
the 1.49 GHz luminosity llCondon et al.l Il99(llf apply¬ 
ing the relatio n betw een Ti .49 GHz and Ti. 4 GHz given by 
ICondon et al . ill 99llf . The obtained value, SFRi. 4 ghz = 
6.4 M 0 yr -1 , is a little lower than those derived from 
previous estimates but, as it was commented before, this 
rate corresponds to a larger time period than the Ha or 
FIR SFR indicators. Using the red dening-corrected Ha 
flux and the expression given bv lDooita et alJ 1 2002 1. we 
can also derive the thermal flux at 1.4 GHz for the galaxy, 
.F 1.4 ghz thermal = 1-94 mj. It corresponds to around 6 % 
of the total 1.4 GHz flux, in agre ement with th e average 
value found in starburst galaxies (Condon 11992). 

Assuming that L r = 6.61 x I Q 10 L m and using 
the relation given by [Gallagher(4984) we derive a 
SFR s = 1.9M 0 vr-U lGalzettil (1200111 give a correlation 
between the SFR of the starburst and the B luminosity of 
the host galaxy. Using the structural parameters derived 
for IRAS 08339+6517 (see TableEJ), we find that the total 
R-luminosity of the underlying component is ~ 15% of 
the to tal. Lb,uc = 8-95 x 10 9 L 0 . Applying the ICalzettil 
l200ll) correlation, we derive SFRs = 1.1 Mq yr _1 . The 
SFR derived from the blue luminosity corresponds to the 
past few billion years, whereas those derived from the Ha 
or FIR fluxes indicate the current (< 10 7 yr) SFR. 

We can also use the relation given bv lKennicutt (jl998h 
between the SFR and the luminosity of [0 11 ]. Using the 
Ha luminosity derived from our images for the entire 
galaxy (see Table n and the [Oil]/Ha ratio from our 
spectra ([Oil]/Ha = 1.01 for the 12" aperture), we find 
SFR[ on] = 13 Mq yr -1 , considering the low-limit of 


the llvennicutf] (1998) calibration, corresponding to blue 
emission-line galaxies. A more realistic relation between 
the SFR and the luminosity of [0 11 ], that take into ac¬ 
count the oxygen abundance and io nization param eter of 
the ionized gas, was presented by iKewlev et al.l ( 2004 1 
(their eq. 15). Applying this relation, we obtain SFR<o n 1 
= 8.9 Mq yr -1 . 

Global soft X-ray luminosity may seem to serve as 
an indicator for the SFR in star-forming galaxies. In the 
last years, several authors have tried to find a correla¬ 
tion between the s oft X-ray emission and the SFR. Using 
the iRanalli et al ] ll2003ll ex pression and the Lx valu e 
from the ROSAT Satellite dStevens fc Strickland! Il998 1. 
Lo. 2 - 2.0 kev = 2.81 x 10 41 erg s _1 , we derived SFRl x = 
61.8 Mq yr _1 , a very high value compare d with that 
found using the others calibrations. IStevens fe Stricklandl 
lll998h find that the X-ray luminosities of WR galaxies are 
considerably higher than those of the other galaxies with 
the same B luminosities, consequence of the higher occur¬ 
rence of superbubbles in WR galaxies. Superbubbles are 
hollow cavities with size of the order of kpc, expanding 
with velocities between 25 and 150 km s" 1 produced by 
the combi nation action o f supernova explosions and stel¬ 
lar winds llLeithereill994lh We can es timate the supernova 
rate, vsn, using the expression (7) in lKewlev et alJ 12000) 
that assumes a Salpeter IMF with lower and upper mass 
limits of 0.1 and 100 Mq and a minimum initial mass for 
supernova detonation of 8 M 0 : 


v SN = 7.42 x 1CT 3 x SFR, 


( 4 ) 


Assuming a SFR = 9.5 Mq yr" 1 , we find a supernova 
rate of vsn ~ 0.07 yr -1 , i.e., a supernova explosion each 
14 years, a bit higher than the va lue usually foun d in 
starburst galaxies [around 0.02 yr -1 dKewle v*et al.l200ci) ]. 
Thus, the relation given bv lRanalli et al.l d2003ll between 
SFR and the X-ray luminosity seems to be not suitable 
for very young starbursts. 

The gas depletion timescale is defined as: 

Tgas = 1.32 x Mhi / SFR, 


( 5 ) 
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Wavelength (A) 


Fig. 15. Part of the spectrum of IRAS 08339+6517 between 
4300 and 5100 A showing the Hell A4686 emission line associ¬ 
ated with Wolf-Rayet stars. The dotted line indicates a linear 
fit to the continuum. 

llSkillman et all2003th indicating the number of years that 
a galaxy may continue to form stars at its current rate. 
The factor 1.32 was introduced to account for He. For 
IRAS 08339+6517 we obtain T gas = 0.15 Gyr (0.68 Gyr 
considering that the mass of the HI tidal tail). It is a very 
low value, indicating the starbursting nature of the galaxy. 

In conclusion, the SFR derived using the Ha, [Oil], 
FIR, 15/itm, 60/irn and 1.4GHz luminosities are in a very 
good agreement and give a SFR ~ 9.5 M© yr” 1 . The 
value of the SFR found using the R-luminosity of the un¬ 
derlying population suggests that the star formation in 
the last 100 Myr has not been as high as the present rate. 
Furthermore, the high SFR derived using X-ray luminosi¬ 
ties implies a strong massive star formation in the last few 
Myr, supporting a high rate of supernova explosions and 
the possible detection of the WR stars. 


these authors ruled this out because they do not detect 
the He n A4686 in their optical spectrum. 

Although the He n A4686 emission line is weak, we have 
used the evolutionary synthesi s models for O and WRpop- 
ulations in young starburst of lSchaerer fc Vacca t l99a) to 
perform a tentative estimation of the WR/(WR+0) ra¬ 
tio. Assuming that all the contribution of the He n A4686 
emission line comes from WNL stars and considering a lu¬ 
mino sity of £ (WNL 4686) = 1.7xl0 36 erg s ” 1 for a WNL 
star llVacca fc Conti 1992il . we find around 310 WNL stars 
in this burst. To derive the WR/(WR+0) ratio, the con¬ 
tribution of the WR stars to the total ionizing flux must 
be considered to obtain the total number of O stars. 
Assuming a l uminosity of L(Hl3 ) = 4.76xl0 36 erg s ” 1 for 
a 07V star llVacca fc Gontilll992lf and n = Q7V/0 = 0.25 
for an age of around 4.5 Mvr Is'chaerer & Vaccalll99^1 . we 
derive around 10700 O stars. This implies a WR/(WR+O) 
ratio of 0.03 in knot #1. An identical value is found using 
the calibration between the W R/(WR+Q) ratio and the 
flux of the WR bump given by ISchaerer fc Vaccal (1998) 
(their equation 17). 

The weak He II A4686 emission line is only detected in 
the spectra of knot ffl (1" aperture) and in the 2.5" aper¬ 
ture but it is not found in any of the other spectra of the 
same galaxy. We suspect that maybe knot #1, which hosts 
the youngest and most powerful burst in the system, is the 
only area where a substantial population may be present. 
Thus, aperture effects and the position of the slit can play 
an important role in the detection 
autho rs have previou sly poin ted 
iLopez-Sanchez e~ Tn i20Q4albh . 

However, although the He II A4686 emission line seems 
to really exist in knot # 1 , deeper spectroscopical data will 
be necessary to confirm that it is broad and/or that the 
WR bump undoubtedly exists. Only under these facts the 
classification of IRAS 08339+6517 as Wolf-Rayet galaxy 
can be established beyond any doubt. 


W R fea tur es, as some 
rng et aTIllc 


4.3. WR population 

IRAS 08339+6517 was included in the X-ray stu dy 
of W R galaxies performed by IStevens fc Stricklandl 
ill 998lf b ecause prev i ous authors llContil 1199 It 
iGonzalez-Delgado et al . Il998h suggested that it has 
recently passed through the WR galaxy phase. However, 
our deep optical spectrum reveals that WR stars could 
still be present in the starburst. In Figure El we show 
part of the spectrum of the central object of IRAS 
08339+6517 (knot # 1 ) between 4200 and 5100 A. It 
shows a very weak bump between 4660 and 4700 A and 
an emission line at A4685.4 that could correspond to the 
Hen A4686 emission line. If this assumption is correct 
it would indicate the first detection of WR stars in this 
starburst gal axy. The UV spectrum of IRAS 08339+6517 
presented by iMargon et all (jl988ll showed a very weak 
emission line around 1640 A that could be attributed to 
Hen A1640, line also associated with WR stars. However, 


4.4. The nature of IRAS 08339+6517 


4.4.1. Radial dependence of metallicity and reddening 

We have analyzed the radial dependence of the de¬ 
rived metallicity of the eight apertures traced in IRAS 
08339+6517 (four apertures corresponding with knots and 
four apertures with different sizes and centered on the nu¬ 
cleus, see Figure^, which is shown in Figurc ITTil The oxy¬ 
gen abu nd ance obtained from the empirical calibration of 
IPilvuginl d200lh (see Tables O and 0 was assumed for each 
aperture. 

Except for knot B, the oxygen abundances derived 
from the nucleus [knot #1 with 12 + log(0/H) = 8.55] to 
the outer regions [knot A with 12 + log(0/H) = 8.42] is 
continuously decreasing. Thus, despite the uncertainties, a 
possible presence of a weak abundance gradient along the 
disk of IRAS 08339+6517 (around —0.043 dex kpc” 1 ) is 
detected. The avera ge gr adient i n spira l galaxies is around 
—0.06 dex kpc ” 1 ( Za/ritski et IDEm 
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Fig. 16. Relation between metallicity and distance from the 
center of IRAS 08339+6517 for the observed knots (triangles) 
and different aperture sizes (open circles). We have consid¬ 
ered the oxyge n abundances derived using the calibration by 
IPilvuginl ( 2001 1. The dotted line is a linear fit to all the data 
except knot B, which shows the higher O/H value. The linear 
ht suggests an abundance gradient along the disk. 


Fig. 17. Relation between the reddening coefficient, C(H/3), 
with distance from the center of IRAS 08339+6517. Observed 
knots are plotted with triangles and the different aperture sizes 
with open circles connected with a dotted line. 

Table 8. Principal properties derived for the ionized gas 
in IRAS 08339+6517 using different aperture sizes. 


Besides the oxygen metallicity, the W a b s calculated 
from the Balmer decrement increases with the aperture 
size, confirming again the existence of an evolved un¬ 
derlying stellar population. The C(H/3) derived for dif¬ 
ferent apertures inside IRAS 08339+6517 also indicate 
an apparent slight gradient between the nucleus and the 
external zones of the galaxy (se e Fi g ure EZL This may 
explai n the differences found by iGonzalez-Delgado et al.l 
: 199S • in the color excess between two different aperture 
sizes. For the nucleus, they estimated E(B — F)=0.19 
(from the spectrum obtained with a 1.7" X 1.7" aper¬ 
ture using the Goddard High-Resolution Spectrograph 
on the Hubble Space Telescope ), similar to the value we 
derive for our 1" x 1" aperture, E(B — F)=0.20±0.01 
[C(H/J)=0.30±0.02]. For the external zones, they derived 
E(B — F)=0.10 (using the 12" circular aperture with 
the Hopkins Ultraviolet Telescope ), which corresponds 
with the one we obtained for our 12" x 1" aperture, 
E{B - F)=0.12+0.01 [C(H/?)=0.22±0.02]. The effect of 
the aperture size on derived inte grated p ro perties of galax- 
ies has been recently studied bv lKewlev et al. 1 ll2005l) . who 
concluded that the difference between the nuclear and 
global derived values for metallicity, extinction and SFR 
is substantial: for example, for a late-type spiral galaxy, 
the metallicity of the nucleus is ~0.14 dex greater than 
the global metallicity of the galaxy. 

4.4.2. Knot B: a TDG or the remnant of a merger? 

As it can be observed in Figure the external knot 
B, located 3.4 kpc from the center of IRAS 08339+6517, 
shows a high metallicity [12 + log(0/H) = 8.58], even a 
bit higher than that derived for the nucleus of the galaxy. 
Note that knots A and B are at rather similar distances 


Aperture (") C(H/3) W nt „(Hi) (A) 12+Iog(Q/H) 


1 

X 

1 

0.30 

± 

0.02 

1.1 

± 

0.1 

8.55 

± 

0.10 

2.5 

X 

1 

0.26 

± 

0.02 

1.4 

± 

0.1 

8.50 

± 

0.10 

5 

X 

1 

0.24 

± 

0.02 

1.5 

± 

0.1 

8.47 

± 

0.10 

9 

X 

1 

0.21 

± 

0.02 

1.6 

± 

0.1 

8.45 

± 

0.10 

12 

X 

1 

0.22 

± 

0.02 

1.8 

± 

0.1 

8.45 

± 

0.10 


from the nucleus but are located on opposite sides of the 
galaxy (see Figure 0 . Their spectra also show approxi¬ 
mately the same S/N. Knot A seems to be a bright Hii 
region (or a complex of H II regions) in the outer areas 
of the disk. Furthermore, as we can see in §3.4.3, knot B 
also shows a decoupled kinematics with respect to the mo¬ 
tion of ionized gas associated to the tidal tail, so it might 
not be symply an intense star-formation region in the out¬ 
skirts of the disk. We consider two possible explanations 
to these facts: 


1. TDG nature. Dwarf objects showing high metal- 
licities and decoupled kinematics are cha racteristics 
of candidate tidal dwarf galaxies, TDGs ( Due et ail 
2!ium. The R-magnitude derived for this region is 
ms ~ 20.1, implying an absolute B -magnitude 
of Mb ~ —14.4. Thus, considering its oxygen 

abundanc e, this object is away from the relation 
given by iRicher fc McGalll ( 11995 ) for dwarf irreg- 
ular galaxies. Fo l lowin g the analysis performed by 
IWeilbacher et al.l (2003), it means that knot B (per¬ 
haps, even all the NW arm, knot #7 in Figure 0 
could be a TDG c andidate. The mean me tallicity of 
the objects in the IWeilbacher et all :V 2()P3 i sample is 
12+log(O/H)=8.34±0.14, but at least five have oxy¬ 
gen abundances higher than 8.55 (see their Figure 3), 
similar to the one derived for knot B. However, a gen- 
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uine TDG must be a self-gravitating entity. We can 
not derive it from our spectrum, so it is not possible 
to confirm that B is really a TDG candidate made 
from material stripped from the internal regions of the 
galaxy. 

2. Merger nature. Knot B could be a remnant of a pre¬ 
vious merging process suffered by IRAS 08339+6517 
in the past. It would explain the disturbed morphol¬ 
ogy in its outer regions, specially the long arc of mate¬ 
rial connecting the north of the galaxy with the south¬ 
ern bright ray (see Figure H) ■ This scen ario was pre¬ 
viously suggested by ICannon et al. i200l and would 
also account for the intensive star-formation activity 
throughout the galaxy. If this early merger really hap¬ 
pened, it would be in a very advanced state. Following 
the detailed analysis o f a galactic me rger sequence pre¬ 
sented bv lHibbard fo Gorkoml ( 1996 1. late-stage merg¬ 
ers have tidal appendages emanating from a single nu¬ 
cleus surrounded by a mostly relaxed stellar profile. 
Their deep R image of NGC 3921, prototype of this 
stage of merging, has morphological signatures similar 
to those found in IRAS 08339+6517. Other interest¬ 
ing minor merger is t he A toms-for-Peace galaxy, NGC 
7252: iHibbard et alJ dl994l) found that most of its Hi 
mass resides outside the merger remnant, also simi¬ 
lar to the case of IRAS 08339+6517. However, if the 
merger hypothesis is right, the interaction with the 
companion galaxy would not be necessary and both 
the optical plume and the detected Hi tidal tail (see 
Figure 113 would be aligned to the direction of the 
companion galaxy only by chance. The effects of the 
interaction with the external companion galaxy can 
not be discarded. Although a minor merger occurred 
at early times the interaction with the external object 
seems to be more important at the present time. 


4.4.3. The non-AGN nature of IRAS 08339+6517 

IRAS 08339+6517 does not host an active galactic nucleus 

(AGN). Several facts support this statement: 

1 . its spectrum, which is similar to the ones of typical 
starburst galaxies; 

2. the FWHM of the emission lines. For example, the 
FWHM of H/3 and [O ill] A5007 are, corrected for in¬ 
strumental broadening, 263 km s -1 and 259 km s -1 , 
respectively; the typical FWHM median value r ange 
for AGNs is between 350 and 550 km s -1 dVeron et a,l.l 
119971) : 

3. the position of t he observational data over the 
IPooita et al . (l2000fl diagnostic diagrams, that are con¬ 
sistent with the loci of typical H n regions and not with 
AGN; 

4. its logarithmic ratio of FIR to radio flux density, q 
= 2.34, that is consistent wit h the valu es derived for 
normal galaxies, q = 2.3 dCondon et alJfl992 i'I: 

5. and the correlation between FIR and radio emission 
that satisfied the galaxy, both using the FIR and 1.49 


GHz luminosities JCondon et alJHl99lll or the 1 .4 GHz 
radio continuum and 60 urn FIR fluxes lYun et alJ 

iS- 


4.4.4. A luminous compact blue galaxy 


The properties observed in IRAS 08339+6517 suggest 
classifying it as a luminous compact blue galaxy (LCBG). 
LCBGs are ~ L* (L* = 1.0 x 1O 1O R 0 ), M B < -18.5, blue 
(R — V < 0.6), high mean surface brightness within the 
half-light radius (SB b < 21 mag arcsec -2 ), vigorous star- 
burst ing galaxies wi th an underlying older stellar popula¬ 
tion dGuzman et al . ljU). Considering the B magnitude 
derived for IRAS 08339+6517, we found that it satisfies 
all these characteristics: it has (B — V) = 0.05 ± 0.08, M B 
= —21.58 ± 0.04 and for 2.6", its half-light radius in the 
R-band, it possesses a mean surface brightness of SB b ~ 
17.1 mag arcsec -2 . 


LCBGs ar e no t common at low redshifts 
( Guzman et al . EhS) and the ir evolutio n and na¬ 
ture is stil l dis cusse d: whereas iKoo et al.1 (1225) and 
iGuzman et alJ (^996) suggested that LCBGs are the 
progenitors of local low-mass spheroidal or i rregular sys¬ 
tems e xperi encing a strong starbu rst phase. IPhillins et alJ 
(ll 997h and iHammer et al.l (l200lh pointed out that they 
are the progenitors of present-day bulges of massive 
spirals. In the local Universe, LCBGs seem to be more 
evolved objects and even they could b e the eq ui valen t 
of the high z Lyman-break galaxies (Erb et al . m, 
being ideal templates for studies of galaxy evolution 
and formation. A recent s ample of LC BGs in the local 
Universe was presented bv lWerk et alJ (2221) who found 
that, on average, LCBGs show strong star-formation 
activity, emit detectable radio continuum flux and have 
lower metal abundances that those expected from the 
luminosity-metallicity r elation fo r starbu rsti n g emi ssion- 
line galaxies given by iMelbourne fc Salzeil ( 200 21. All 
these properties are also satisf ied by IRAS 08339+6517 
[the IMelbourne fc Salzerl (’ 2002 1 relation gives an oxygen 
abundance of ~9.2, much higher than that we derived, 
12+log(0/H) = 8.45+0.10], 


4.4.5. The nature of the dwarf companion galaxy 

The dwarf compa n ion ga laxy follows the relation given by 
iRicher fc McCali ( Il995fl between absolute R-magnitude 
and oxygen abundance for dwarf irregular galaxies, indi¬ 
cating that it is not a TDG formed from material striped 
from IRAS 08339+6517 but an external and independent 
dwarf galaxy. The apparent component of solid-body ro¬ 
tation found in Hi, the indications that the Hi gas of 
the tidal tail has been probably stripped from the main 
galaxy and the lower star-formation activity found in the 
companion galaxy suggest that the effects of the interac¬ 
tion are not as intense as in IRAS 08339+6517. 
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Fig. 18. Deep R image of IRAS 08339+6517 and its compa nion g alaxy showing the faintest details. Our continuum-subtracted 
Ha image (white contours) and the Hi map obtained by lCannon et ah (1200 41 (black contours) is superposed. Note that the 
weak optical plume in the direction of the companion galaxy coincides with the HI tidal tail. 


5. Conclusions 

We have used deep optical and Ha imagery together deep 
optical intermediate-resolution spectroscopy to analyze 
the morphology, colors, ages, stellar populations, physi¬ 
cal conditions, kinematics and chemical abundances of the 
galaxy IRAS 08339+6517 and its dwarf companion galaxy, 
2MASX J08380769+6508579. Our data r einforce the re- 
sults of the HI observations performed by ICannon et a,l. I 
J2004 that both objects are in interaction and that the HI 
tidal tail seems to be formed mainly by material stripped 
from IRAS 08339+6517. 

We have obtained the oxygen abundances of both 
galaxies using empirical calibrations. The O/H and N/O 
ratios of the two galaxies are rather similar, suggesting 
a similar degree of chemical evolution of the system. We 
have analyzed the chemistry of the ionized gas using dif¬ 
ferent aperture sizes and found an apparent metallicity 
gradient between the nucleus and the external zones of 
IRAS 08339+6517. The reddening coefficient also shows 
differences between the central and external zones. 

IRAS 08339+6517 shows important Ha emission in its 
inner regions. IRAS 08339+6517 is not powered by an 
AGN but it is a nuclear starbursting galaxy. We have de¬ 
tected at least two different populations in the galaxy, the 
age of the youngest one being around 4 6 Myr. But 

the more evolved stellar population, with age older than 
100 - 200 Myr, fits an exponential intensity profile. A 
model which combines 85% of a young (6 Myr) popula¬ 
tion with 15% of an old (140 - 200 Myr) population can 
explain both the spectral energy distribution and the H 
Balrner and Hei absorption lines observed in our spec¬ 
trum. Furthermore, the NIR colors suggest that even an 


older stellar population, with age upper than 1-2 Gyr, is 
also present in IRAS 08339+6517. Its companion galaxy 
is practically dominated by an old population with age > 
250 Myr, although it also hosts a recent ~ 6 Myr starburst 
in its external areas. The Keplerian, ionized and cold dust 
masses derived for the main galaxy seem to be the usual 
for young starbursts. 

Our spectra seem to show weak Wolf-Rayet features in 
IRAS 08339+6517, but they are located only in the central 
knot of the galaxy. However, deeper spectra are needed to 
confirm the WR nature of this galaxy. We discuss that 
aperture effects and localization of the bursts with WR 
stars could play a fundamental role in the detection of 
this sort of massive stars in starburst galaxies. 

We have derived the SFR of IRAS 08339+6517 using 
multi-wavelength correlations. The Ha, FIR and radio lu¬ 
minosities give similar values, around 9.5 M 0 yr , sug¬ 
gesting that the contribution of reddening by dust has 
been properly estimated. The high SFR derived using X- 
ray luminosities implies massive stars formation in the last 
few Myr, as seems to be supported by the high rate of su¬ 
pernova explosions and the probable detection of a large 
number WR stars. 

Knot B, located in the NW arm of the main galaxy, 
has the higher metallicity and possesses peculiar kinemat¬ 
ics. It could be a TDG candidate because its metallic¬ 
ity is higher than that expected for a dwarf galaxy with 
its R-luminosity and because the kinematics of the ion¬ 
ized gas are decoupled from the general kinematic pattern. 
However, it could also be a remnant of a previous early 
merger because of the peculiar morphology of the galaxy 
and the existence of the Hi tidal tail. Besides it, at the 
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present time the interaction with the external companion 
galaxy is much more prominent and probably the origin 
of the HI tidal tail. 

Finally, IRAS 08339+6517 could be classified as 
a luminous compact blue galaxy (LCBG) because of 
its color, absolute magnitude and surface brightness. 
There are very few local LCBGs nowadays detected but 
nearly half of them have optical companions, present 
disturbed morphologies and/or are clearly interacting 
1 Garland e t~m 12QQ-4 An example o f local LCBG is 
Mkn 1087 llLopez-Sanchez et al.l2004bl) , that is in interac¬ 
tion with two nearby galaxies and shows bridges and tails 
that connect dwarf surrounding objects with the main 
body. If interactions were the responsible of the activity 
in LCBGs, it would indicate that they were perhaps more 
common at high redshifts, a s the hierarchical galax ies for¬ 
mation mod els predict [i.e., iKauffmann fc White! ( 1993 0: 
ISoringel et al.l |2005)]. This fact would also support the 
idea that interaction with dwarf companion objects could 
be an important trigger mechanism of the star formation 
activity in local starbursts. 
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